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The instantaneous distribution of chain compositions and chain lengths in vinyl copolymers 
is obtained in a simple form valid for long chains. The compositions of chains of a given length 
are normally distributed about the mean value, with a standard deviation which can be calcu- 
lated from experimentally observable quantities. The distribution of chain lengths intimately 


resembles that for simple polymers. 





I. INTRODUCTION 


ECENTLY several authors'~* have inde- 
pendently derived a simple differential 
equation! relating the over-all composition of a 
long-chain vinyl copolymer to that of its parent 
monomer solution, and this relation has been 
found to agree well with experimental data.? 
The distribution of sequences of like monomer 
units along the chains was also readily ob- 
tained.! 3 
Since the molecular weight of a copolymer is 
finite, the compositions, as well as the lengths, 
of the individual molecules cannot all be iden- 
tical; even in the copolymer formed during a 
very small time interval, there exists a distribu- 
tion of chain lengths and compositions. Although 
one could reasonably expect the fluctuations of 
composition to be quite small for long molecules, 
it seems worth while to obtain a quantitative 


* Contribution No. 513 from the Research Labora- 
tory of Physical Chemistry, Massachusetts Institute of 
Technology. 

1 T. Alfrey and G. Goldfinger , J. Chem. Phys. 12, 205 and 
322 (1944). 

2F. R. Mayo and F. M. Lewis, J. Am. Chem. Soc. 66, 
1594 (1944). 

3F. T. Wall, J. Am. Chem. Soc. 66, 2050 (1944). 

4 See Eq. (5). 


estimate of their magnitude. A general theory of 
copolymerization which embraces this problem 
has been formulated by Simha and Branson,* but 
the equations given by them are unfortunately 
too complex for convenient use.® With the aid of 
further approximations (but no new assump- 
tions) valid for long chains, the distribution of 
chain compositions and lengths can be expressed 
in a simple and tractable form. In this paper, 
the simplified result is derived, discussed, and 
illustrated. 


II. MEAN COMPOSITION AND CHAIN LENGTH 


In describing the mechanistic assumptions and 
the symbolism to be employed, some repetition 
of the discussion in previous papers!~* inevitably 
occurs. We consider a binary system with 
monomer concentrations A and B. The growing 
free-radical’? chains present during polymeriza- 


5R. Simha and H. Branson, J. Chem. Phys. 12, 253 
(1944). 

6 Some of the equations and all of the numerical examples 
in the paper of Simha and Branson pertain to a special case 
which rarely, if ever, is found experimentally. 

7 The expression ‘free radical’ is used advisedly, in 
conformity with present ideas regarding the most common 
polymerization mechanism. Actually the equations to be 
derived are valid for any mechanism in which the concen- 
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Fic. 1. Reduced composition distribution according to 
Eq. (17). Weight fractions are plotted against the reduced 
composition-deviation at relative lengths of 0.25, 1, and 4. 


tion are called either A-radicals or B-radicals, 
according to the identity of the monomer unit at 
the active end. The total concentrations of A- 
and B-radicals of all lengths and compositions 
will be written m and n, respectively, while the 
symbols m,, and n,;, denote the concentrations of 
such radicals composed specifically of r A-units 
and s B-units. If, as is usually assumed, the 
reactivity of a radical depends only on the 
nature of its active end (and not on its length or 
composition), just four different propagation 
(growth) reactions must be considered,® as set 
forth in,the following scheme: 


(RaaM,sA), 
Mre+B—ny,s41 (Ravm,sB), 
Nr+A—mrs1,2 (RoatreA), 
(Rutt B ). 


Mrz t+A—M,41, 8 
(1) 


Nr+B—n,, s+1 


The expressions in parentheses following each 
equation denote the reaction rates. The bimo- 
lecular character of these propagation steps may 
be regarded as experimentally confirmed,? and 
no other specific kinetic assumptions are required. 

The rates of formation of the two kinds of 
radicals by other than propagation steps will be 
denoted simply by 7,4 and %B. These rates thus 
tration of active growing chains (whether free radicals, 
ions, or excited molecules) is limited by some termination 
reaction to such small values that the steady-state as- 
sumption is valid. 

8 the aforementioned articles, references 1-5, or 
earlier papers by R. G. W. Norrish and E. F. Brookman, 


Proc. Roy. Soc. (London) A171, 147 (1939) and E. Jenckel, 
Zeits. f. physik. Chemie A190, 24 (1941). 
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Fic. 2. Reduced chain-length distribution according to 
Eq. (17). Weight fractions are plotted against relative 
length at reduced composition-deviations of 0, 0.5, 1, and 
1.5. The dashed curve is the distribution summed over all 
compositions, divided by two to facilitate comparison. 


include not only true initiation processes, but 
also the possible formation of A- and B-radicals 
by addition of monomers to other types of 
radicals produced in chain transfer. Obviously, 
4, and % are therefore usually not constants, but 
are functions of monomer, solvent, and catalyst 
concentrations. Similarly, the rates at which the 
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Fic. 3. Over-all distribution of compositions, according 
to Eq. (14), for copolymers with \/20g0= 400, showing the 
effect of the parameter x. 


® For example, for the simple case of a peroxide-activated 
polymerization in which chain transfer is negligible, one 
might expect the following expressions to hold: 1¢4A =kiaCA, 
4B=kpCB, ram=Rksmm?+k32mn, ren=ksrzmn+k3nn’, in 
which C is the peroxide concentration and the various k’s 
are rate constants. 














CHAIN LENGTHS 


two kinds of radicals are converted into stable 
polymer (including chain-transfer processes) will 
be written 74m and 7), where the quantities 
7, and 7 are certainly not constants.°® 

With these definitions, the total radical con- 
centrations are governed by the equations 


dm/dt=1,A —kawmB+kyanA —tam, 
dn/dt = B = kpanA ob kamB — Ton. 


In practical cases, the chains are so long that 
propagation must be much more rapid than 
initiation or termination. With the neglect of 
these latter, the assumption of a steady state! 
(dm/dt=dn/dt=0) then gives the important ap- 
proximate equality 
kamB = kpanA P (3) 

which, together with the rates of monomer con- 
sumption, 

—dA/dt=RaamA+kyanA, 

—dB/dt=kynB+kanB, 
leads to the previously mentioned expression!~* 


for the over-all composition of the copolymer 
which is forming at any chosen time: 


dA/dB=A(aA+B)/B(A+£B) =po/qo, 


where 


(2) 


(4) 


(5) 


a=Raa/Rka, B=kw/Roa- 


For later use, we have here introduced the 
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symbol pp>=1—go for the over-all mole fraction 
of A in the copolymer chains produced during 
the (infinitesimal) time interval when the mono- 
mer concentrations are A and B. The convenience 
of dealing with symmetrical equations leads us 
to retain a separate symbol go for the over-all 
mole fraction of B. The ‘“‘monomer reactivity 
ratios’? a and 8, which with the aid of Eq. (5) 
may be obtained from experimental studies of 
copolymer compositions, also play a role in de- 
termining the distribution of chain compositions. 

The number-average length (polymerization 
degree) X of the growing radical chains is 
simply the total rate of monomer consumption, 
—d(A+B)/dt, divided by the rate either of 
radical production or of radical termination ;" thus, 


—d(A+B)dt aA?+2AB+ 6B? 
ram+rn (tA ka) + (7oB/ Roa) 


with the aid of Eqs. (3) and (4). If all the chains 
terminate by disproportionation of radicals, the 
number-average length of the stable polymer 
chains is also ); if all the chains combine in pairs 
when mutually terminating, it is 2A. In view of 
the unspecified concentration-dependence of the 
termination parameters tr, and 7, Eq. (6) is of 
no direct use as it stands. We shall employ it 
only to express the chain-length distribution con- 
veniently in terms of the mean chain length. 





III. DISTRIBUTION OF COMPOSITIONS AND LENGTHS 


When the steady-state condition is applied separately to radicals of each length and composition, 


the following recursion formulae result: 


dm,. 
_ RaaA M,y—1, ot RiaA Nr—1, 
dt 


anys 


= kwBn,, si tkwBm,, s—1 


ae (RaaA +kaB +1a)Mrs = 0, 
(7) 
— (kwB+ kiaA +70) trs=0. 


The species m9 and 1, of course, follow different equations, involving i,A and iB, respectively. 


The solution of Eqs. (7) is® 


r é. r—1 s—1 


Mrs = Wa" wr* 
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j-\ j 


sae (nen) brief discussion of the conditions of validity of this assumption, see W. H. Stockmayer, J. Chem. Phys. 12, 
11 On the steady-state assumption, radicals are made and destroyed at equal rates, so that i,A +iB=ram+ron. 
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where 


Wa= RaaA /(RaaA +kwBt+ra), b= kwB/ (k»B +RiaA +7»), x= Ravkba/Raakoe = 1/ap. 


For the pure radical species m,o and mo, the equations are simpler, but they need not be considered 
here. 

A physical interpretation of Eqs. (8) is not difficult.’ Each expression falls into two sums, the 
first describing chains which began with an A-unit and the second those which began with a B-unit. 
A single term within a sum denotes radicals in which 7 cross-additions of one kind (A to B in the first 
equation and B to A in the second) have occurred, and the combinatory factor gives the number of 
ways (sequences of addition) in which such radicals can be built up. The fractions w, and w are 
propagation probabilities for homogeneous addition (A to A and B to B), while x is a measure of 
the relative probability of cross-addition. 

The procedure followed in simplifying Eqs. (8) is described in some detail in the Appendix. In 
brief, the sums over j are approximated by integrals and evaluated with the aid of Stirling’s formula. 
The logarithm of m,, can then be expanded as a power series in the composition-deviation y, de- 
fined by 

y=p—Po=q-4Y, (9) 


where p=1—g=r/(r+s) is the composition of an individual chain and pp>=1—gp is the over-all 
composition. This procedure leads naturally to a Gaussian distribution” of chain compositions about 
the mean value. To this approximation, the distributions for A- and B-radicals are identical. 

The final result is most compactly expressed as a continuous weight distribution, in which equal 
weights are assigned to the two types of monomer unit. The weight fraction of radicals with lengths 
(polymerization degrees) between / and /+d/ and composition-deviations between y and y+dy is 
then given by 





aioe w(I, y)dldy=[exp (—1/d)ldl/d* JL (l/2epogox)* exp (—Ly?/2pogox)dy ], (10) 
k=[1—4pogo(1 — af) }? 
aA?+2a8AB+ BB? 
~ @A?4+24B+6B? | 


and d is the number-average length of the radicals given in Eq. (6). 

If the termination step is entirely one of disproportionation, Eq. (10) also describes the distribution 
of the stable polymer forming at any instant, for then the concentration N,, of stable chains of the 
indicated length and composition is determined by the equation 


dN,s/dt=TaMrst+ToMrs. (11) 


More generally, if a fraction p of the radical chains were terminated by combination, the stable 
polymer distribution would be 


w(I, y)dldy=[exp (—1/)(1—p+pl/2n)ldl/d* | XL (l/2mpogox)! exp (—Ly*/2pogox)dy], (12) 


in which X still denotes the number-average length of the radicals. This result is in complete analogy 
with the size distribution for pure polymers.'* The radical distribution, of course, is unaffected by 
the nature of the termination process, being always given by Eq. (10). 

It should be emphasized that the distribution equations are as generally applicable as the experi- 
mentally confirmed Eq. (5), for no additional assumptions are required for their derivation. Only 


2 For copolymers with a very small percentage of one component (pogo), an equivalent Poisson distribution could be 
obtained if desired. 
13 E. F. G. Herington and A. Robertson, Trans. Faraday Soc. 38, 490 (1942). 
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the uncertainty with respect to the termination process mars their quantitative utility, and this can 
be seen to affect only the distribution of lengths. It may be concluded that Eqs. (10) and (12) are 
no less valid than the theoretical distributions heretofore derived for pure polymers. In the para- 
graphs to follow, we shall for simplicity use Eq. (10) as it stands. 

To find the chain-length distribution irrespective of composition, Eq. (10) must be integrated 
over y. The limits of integration may be taken as infinite because the distribution of compositions is 
narrow, and in fact the equation is already normalized to these limits. The result is therefore simply 


(w(1))dl=exp (—1/d)ldl/d?, (13) 


which is exactly the equation for pure polymers terminating entirely by disproportionation."® 
The over-all distribution of compositions irrespective of chain length is found in similar fashion by 
integrating over /, the result being 


3dz 
(w(y))dy =———_—, (14) 
4(1 +27)5/2 


2?=)y?/2pogox. 


It is interesting to express the parameter x, which controls the spread of chain compositions, in 
somewhat different terms. For sufficiently long chains, the number-average lengths of sequences of 
like monomer units are given’? by 


l4=1+(aA/B), lg=1+(8B/A). (15) 
Employing these relations together with Eq. (5), one obtains 
k=po(ls—1)+qo0(la—1), (16) 


identifying x as a certain kind of average sequence length. Thus longer sequences correspond to a 
broader distribution of compositions, as would be expected. 

For the special ¢ase a8=1, which corresponds to an earlier theory of Wall,!4 and which also was 
given great prominence in the work of Simha and Branson,® Eq. (10) shows that «=1. In fact, for 
this case the distribution of compositions could have been written down directly, for if a8=1 the 
probability of chain growth depends solely on the nature and concentration of the monomers, and 
not at all on the type of active end possessed by the growing chain. 

The distribution may be more compactly expressed in terms of the relative length, L=//\, and 
the reduced composition-deviation, z= (A/2pogox)y, which has already been introduced in Eq. (14). 
With these substitutions, Eq. (10) becomes 


w(l, y)dldy= W(L, 2)dLdz=x-L! exp [—L(1+2") ]dLdz. (17) 


The shape of the W(L, z) surface thus defined can be visualized from Figs. 1 and 2, in which W(L, z) 
is plotted against z for several values of L, and against L for several values of z. It is seen that greater 
deviations of composition occur in the shorter chains, and conversely that chains with the average 
composition are the longest. 

The magnitude of the effect produced by a change in x may be judged from Fig. 3, in which the 
over-all distribution of chain compositions, given by Eq. (14), is plotted for several systems with an 
arbitrarily chosen value of \/20¢o=400. The abscissa scale is seen to be proportional to the quantity 
(2pogox/d)!, which is a convenient measure of the breadth of the distribution, for integration of 
Eq. (14) shows that 88.4 percent by weight of the copolymer falls in the range |y| <(2poqox/A)}. 

It must be emphasized that the practical utility of Eq. (10) and its corollaries is severely limited, 
in that only the distribution for the copolymer forming at any instant is described. Since the 


with 


4F, T. Wall, J. Am. Chem. Soc. 63, 1862 (1941). 
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parameters Po, go, x, and d will in general change as the reaction proceeds,'!* these equations may 
safely be used only for copolymers of very low conversion. For higher conversions, the instantaneous 
spread of compositions will often be quite insignificant in comparison with the steady change of the 
mean composition, governed by Eq. (5), as the reaction progresses. In principle, the complete dis- 
tribution for a high-conversion copolymer could be found by integration from Eqs. (5) and (10) 
and a detailed knowledge of the dependence of \ on concentration and conversion. In the absence of 
the latter information, such an effort would at present necessarily be speculative and scarcely 
worthwhile. 

It should be possible to subject the distribution formulae to experimental test, preferably by 
fractionating ‘‘azeotropic’’ copolymers of low conversion, These should have rather low average 
molecular weights, in order to secure a sufficiently broad distribution of compositions. 

Because of their simplicity of form and interpretation, these distribution equations should be of at 
least qualitative assistance to a better understanding of the relationship between structure and 
properties of copolymers. 


APPENDIX 


Since the symmetry of Eqs. (8) is obvious, it will be sufficient to consider the expression for m,,’ 
which can be rewritten in the form 


Mrs = (wa?ee?) "[(ta/ Raa) (So/g—S1/pq) + (to/ Rav) (S1/ Pq) J, (A1) 
where 
l=r+s, p=r/l, q=s/l=1-p; 
and 
Sn= 5 (G/l)"T;, 
with 


T- asosi(”)( ’). 


For reasons of symmetry, it is convenient to retain both p and gq as symbols. 
For large /, Stirling’s formula may be used to evaluate the combinatory factors. Putting j=/w, 


we have 


and hence 


\(). on 
In T ;2—I1n (221) —3 In [(p—w) (q—w)/bq] 
+I[w In x—2w In w—(p—w) In (b—w) —(q—w) In (g—w) ]. (A3) 


By differentiating this expression, we obtain 


1 1 _ — 
0 In Ty/ow=3( + )+! In = w)a “| (A4) 


p-w q-w w 








from which it is seen that T; passes through a maximum at a certain value of w. For large /, only the 
last term in Eq. (A4) need be retained. Setting wmsx =u, we therefore have 


u?=x(p—u)(q—u), (AS) 


u={1—[1—4pq(1—1/x) ]'}/2(1—1/x), (A6) 


where the negative root is necessary because “(=Wmax=jmax//) may not exceed unity. 


15 In the special case of an “‘azeotropic”’ copolymerization, reference 3, po and « remain constant, but A may still change. 
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We may pause to recall that 7 is the number of cross-additions of one kind (actually, A to B) in 
the chain. Hence u(=jmax//) is the most probable fraction of such bonds in a chain of length /.and 
composition p. For chains with the mean composition po, the use of Eq. (5) of the text, together with 
the relation x=1/a8, gives, from Eq. (A6), the result 


uo=AB/(aA?+2AB+6B?), (A7) 


which is readily shown® to be the fraction of A— to —B additions taken over all the chains. 
The value of a general term 7; may now be found in relation to the maximum term T,, by setting 
w=u-+v and expanding in powers of v. With the aid of Eq. (A5) to effect cancellation, we find 


(1/2) In (7 j/Tm) = — pv? +O(v’, 1/2), (A8) 


with w?=(1/u)+(1—2u)/2(p—u)(q—u). 

The Euler-Maclaurin sum formula shows that, for large /, the sums defining Sy and S; may be 
replaced by integrals with little error. Furthermore, the limits of integration may be extended to 
infinity, because the maximum in 7’; is sharp; in fact, as is usual in such problems, the use of infinite 
limits is consistent with the previous approximations. We therefore may write 


+00 
So= Li T2 {TAIT f exp (—ly?v?)dv= (al)! Tin/u. (A9) 
Since the maximum term in S; occurs at practically the identical value of 7, we also have 
+0 
Si= Li wT IT f (u+v) exp (—ly?v?)dv=uSp. 


At this stage, Eq. (A1) can be rewritten to read 


Mrs = (wa? wp?) sf (ta/Raa) Fy. + (t5/ Ra) (w/a) (Al 1) 
q 


but the dependence of S» on / and p, g must be more explicitly stated. By substitution for 7, and yu 
in Eq (A9), with the aid of Eqs. (A3), (A5), and (A8), there is obtained 


Pq 2 1—2u 
a 
(p—u)(q—u) u (p—u)(q—u) 


p 
+ p In ( )+9 In ( ; )I (A12) 
p-Uu q-—u 
for chains of length / and composition p, q. 


We next proceed to expand In Sp about its value for chains having a composition fo, go which 
equals the over-all composition as given by Eq. (5) of the text. The deviation from this mean com- 
position will be defined as 





In Sy= —3 In (2ml) +3 In 





y=p—Po=Go—4. (A13) 
Before attacking Eq. (A12), we require an expansion for u. Eq. (6) leads to 
u —Uo= Y(go— po) /K+¥"L(qo— Po)*(1—1/x) — x? ]/n?+O(y’), (A14) 
where the important new quantity x, which persists in the final equations, is defined by 


k=[1—4pogo(1—1/x) ]} 
= (aA?+2a8AB+£B?)/(aA?+2AB+ 6B’). 
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Now, upon entering Eq. (A12) with Eqs. (A13) and (A14), there is obtained, after a lengthy but 
straightforward expansion and collection of terms, the result 


2 pl ql 
s= (=) ( = ) ( = ) exp (—Ly?/2pogox), (A16) 
2nlk Po— Uo qo— Uo 


in which terms of O(y*) in In .Sp have been neglected. 
Before inserting this expression into Eq. (A11), we note the equalities 





Po/(po— uo) = (aA +B)/aA = (RaaA +karB)/RaoA, 
(A17) 


qo/(qo— Uo) = (BB+A)/BB= (kw B+ kiaA)/kwB, 


which may be obtained from Eqs. (5) and (A7). Also, the dependence on composition is so over- 
whelmingly dictated by the last factor in Eq. (A16) that everywhere else the arbitrary composition 
bp, g may be replaced by the mean composition fo, go without appreciable error. It can be shown that 
this maneuver amounts to neglecting a skewness in the distribution which becomes negligible for 
long chains. With the aid of the definitions of w, and w, given in Eq. (8) of the text, Eq. (A11) then 
yields 


A(igA +i,B) 
Mrs = 
ka(aA?+2AB+ 6B?) 


: ( RaaA+kaB )( ky Bt RiaA ) 


£'(2mpogoxl)—* exp (—Ly?/2pogox) (A18) 











kaaA+kaB+ta kwB+ kraA +7» 
The significance of the two factors in as propagation probabilities, for A— and B-— radicals re- 
spectively, is clear. 

To obtain the total concentration of radical chains, summation over all compositions and lengths 
is required, but the approximations previously introduced are obviously consistent with integration 
over ly, followed by summation over /. Eq. (A18) is therefore more conveniently used in the form 





A (igA +%B) l 
mi(y)dy= ( 


at 
xp (—Ly?/2 dy, A19 
balaA?+2AB4+6B") ) exp (—Ly?/2pogox)dy (Ai9) 


2mpogok 


where m:(y)dy is the concentration of A-radicals with length / and compositions lying between 
(potty) and (po+y+dy). From the symmetry of Eqs. (8), a similar equation obviously applied to 
the B-radicals, n:(y)dy being given by Eq. (A19) with the factor A/ka replaced by B/kpa. Thus the 
steady-state condition of Eq. (3) is seen to hold not only for the total radical concentrations, but also 
separately for radicals of each length and composition. 

The concentration of radicals of length / but of all compositions is found by integration of Eq. 
(A19) to be 


+00 
_— f m,(y)dy =A (inA +4,B)E"/ qo (aA?-+2AB+6B?). (A20) 


—@ 


Summation over / then gives the total concentration of A-radicals, yielding 


m= >) m,=[(A (icA +%B)/kav(aA?+2AB+ 6B?) \[t/(1—£) ]. 
121 
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For long chains, must be nearly unity, and a sufficient approximation may be found by expanding: 


£=(1+70/(kaaA +hasB) -™1+70/ (kw B+ kvaA ) |-% 
=1—[pota/(RacA +kaB) |—(gors/(kinB+RraA) ]+--: (A22) 
=1—[(taA /Rar) +(16B/Roa) ]/(aA*+2AB+6B*)+---. 
Using this relation in Eq. (A21), we obtain 
m=(A/kav)((taA +%B)/[(taA /Rav) +(1+B/Rra) }), 
n=(B/Rta)((taA +%B)/[(taA /Rar) +(7+B/Roa) }). 


and similarly 


These results are seen to be consistent with the steady-state requirement that 
1A +%B=ram+7on. 
The number-average length of the radicals is at once found to be 
A=D Imi /L m= LV Ini/YL m=1/(1—-&) 
= (aA?+24AB+6B?)/[(1aA /kav) +(1oB/Roa) J; (A26) 


as previously found more directly in Eq. (6) of the text. 

We finally convert Eq. (A19) into a weight distribution. If equal weights are assigned to A- and 
B-monomers, the weight fraction w:(y)dy of chains of length / and compositions between (po+y) 
and (po+y+dy) is given by 


wi(y)dy =lmi(y)dy/>. lm, 
=1(1— &)?'"1(1/2epogox)* exp (—Ly?/2pogox)dy. (A27) 


Since long chains have been assumed, the distribution of chain lengths can also be conveniently 
made continuous. With the substitution \=(1—£)—! from Eq. (A26), the result is 


w(l, y)dldy=[exp (—1/d)ldl/d? J (l/2mpogox)* exp (—Ly?/2pogox)dy ], (A28) 
which is Eq. (10) of the text. 
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Target Chemistry of Zirconium and Yttrium; Isolation of Zr*® in Pure State 


M. H. KursBatov AND J..D. KuRBATOV 
The Ohio State University, Columbus, Ohio 


(Received January 16, 1945) 


Radioactive Zr®® obtained by bombardment of Y.O; with deuterons, was separated without 
addition of stable zirconium in quantities ~10~ gram. Oxychloride solutions of Zr®* mixed 
with inactive Zr were used for coagulation studies in the concentration range 1.5 10-*— 
1.5X10—! gA/l. It was found that the fraction coagulated increased with dilution under the 
constant conditions used (pH—4.0; NaCl—2.5x10-*?N; T—25°), to 2.7X10-® gA/l. The 
log-log plot of zirconium concentration (C) in solution against coagulated part (c) closely ap- 
proaches a straight line with exponent »=1.50 and k=490. At 1.5X10-® gA/I and lower 
concentrations the coagulated fraction reaches a maximum and becomes constant, k= C/c=.45. 
Within the range studied salt increases coagulation of zirconium above 10-5 gA/I; but below 
10-6 gA/I the salt effect is reversed. Increase in pH results in increase in the fraction coagulated 
(r) in all concentrations studied. At a minute concentration r= (pH —K)/m, K being 1.85 with 
NaCl, 0.05. and 1.58 with NaCl, 0.005N and m being 3.15. The separation of Zr®® in minute 
quantities without addition of stable zirconium is plausible if a Nernst distribution of zirconium 
between coagulated part and solution is postulated. The concept of the solubility of the hydrous 
oxide of zirconium as a constant quantity remaining in solution under specified conditions was 


not verified, 





I, INTRODUCTION 


SOLATION of artificial radioactive isotopes 
without dilution by inactive isotopes of the 
same charge number is pertinent to the solution 
of some problems in nuclear physics. Isotopes of 
a different charge number (Z+1, 2) than those 
bombarded (Z) are produced by several nuclear 
reactions. However, the yields are minute quanti- 
ties so their separation involves unconventional 
procedure and the little known behavior of 
elements in extreme dilution. 

For instance, it has been shown in the past!~® 
that some natural radioactive isotopes which 
form slightly soluble hydrous oxides can be 
agglomerated at very low concentrations even 
when the amounts in solution are much below 
their solubilities. These agglomerates were termed 
radio-colloids while the phenomenon was treated 
mainly as adsorption in order to avoid conflict 
with the solubility product principle.’ 


1 Blau and Rona, Akad. Wiss. Wien 139, 275 (1930). 

2 Chamie, J. Phys. Rad. 10, 44 (1929). 

3 Chamie and Guillot, Comptes rendus 190, 1187 (1930). 

‘Hahn, A — Radiochemistry (Cornell University 
Press, Ithaca, New York, 1936). 

5’ Hahn and Werner, Naturwiss. 17, 961 (1929). 

6 Werner, Zeits. f. physik. Chemie A156, 89 (1931). 

7In the literature there is considerable controversy on 
radio-colloid formation. One school of thought admits only 
an adsorption grouping while another supports independent 
coagulation. An example which has had contradictory 
explanations is the formation of radio-colloids of thorium 
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In the field of artificial radioactivity the for- 
mation of radio-colloids has been used in the 
separation of yttrium isotopes (Y* 8%) obtained 
by the reactions, Sr*® §*87(d, m and d, 2n)Y*% 88, 
without addition of inactive yttrium as carrier.® 
It was found that partial agglomeration of 
yttrium takes place in solutions of minute* con- 
centrations and its extent depends primarily 
upon two factors, the pH of the solution and the 
salt concentration. These two factors act in 
opposite directions and affect the splitting of 
yttrium betwéen the solution and the radio- 
colloid. 

The experiments showed furthermore, that 
with very dilute solutions the fraction coagulated 
at a given pH and salt concentration, increases 
with dilution to some limit. This, of course, is 
contrary to the behavior of yttrium in con- 
ventional concentrations. 

The present work on the tetravalent element 
zirconium was undertaken as an extension of the 
previous studies of yttrium to confirm the basic 


B, an isotope of lead, agglomeration of which has been ob- 
served in concentrations as low as 10-" gA/liter while the 
solubility under the conditions of the experiment was 
5X 10-5 gA/liter. Paneth, Mitt. Rad. Inst. Wien 34, 47, 55 
(1913); Zeits. f. Electrochemie 31, 5 (1925). 

® J. D. Kurbatov and M. H. Kurbatov, J. Phys. Chem. 
46, 441 (1942). 

* Minute quantities are accepted arbitrarily as those less 
than 10-* gram. 
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observation of increased coagulation with dilu- 
tion, within certain limits. The radioactive iso- 
topes of zirconium were selected mainly because 
of the convenient half-lives and strong radiation 
intensities of Zr®® and Zr. 

Of these, Zr® has characteristics which make 
it advantageous for studies of dilute solutions 
(10-* to 10-* gram atom per liter) and Zr*®® for 
concentrations below 10~-° gram atom per liter. 

The interrelation of the isotopes of yttrium, 
zirconium and columbium is shown in Fig. 1. 

The half-lives (J) and upper limits of radiation 
energies in Mev are given according to latest 
data.*!° The radioactive isotopes used in the 
present work are shown in heavy lines. 


Il. TARGET CHEMISTRY OF ZIRCONIUM 


By activation of metallic zirconium with 
deuterons in the Ohio State University cyclotron, 
the radioactive isotope Zr® of T—65.5 days was 
secured for studies of dilute solutions. 

In spite of the convenient half-life and strong 
radiation intensity Zr® has the disadvantage of 
chain disintegration into Cb**, 7— 38.7 days, the 
latter an isomer of stable Cb. Therefore, a 
separation of Cb from Zr is required just prior 
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Fic. 1. Region of yttrium, zirconium, and columbium 
isotopes. Stable isotopes are shown in rectangles and 
radioactive in ovals. 


to each series of experiments with dilute solu- 
tions of zirconium. 

The target chemistry of the metallic zirconium 
(ca. 0.1 g) bombarded with deuterons consisted 
in dissolving the metal, separating the isotopes 
of columbium (d, m and d, 2n reactions) and 


°M. L. Pool and J. E. Edwards, Bull. Am. Phys. Soc. 19, 
7 (1944). 


0G. T. Seaborg, Rev. Mod. Phys. 16, 13 (1944). 


yttrium (d, a reaction) and removing from the 
zirconium the usual active impurities Cu, Ni, 
Co, Na, etc. 

After only the surface of the zirconium was 
dissolved in acids, stable yttrium was added and 





ELECTROMETRIC TRATION 
pH ZrOGi, with 0.01 NNGOH 
60 
- copayiation ——> 
40 
30 

cence 

225 35 «8 05 05 15 125 135 


NOOH g mol. x10® 


Fic. 2. Electrometric titration: final concentration of 
zirconium 4.4 107-3 gA/I. 


later separated by the a acid pro- 
cedure. 

The columbium isotopes were first separated, 
without addition of stable columbium, by means 
of the potassium carbonate procedure. The 
minute quantities of potassium hexacolumbate 
so obtained were used later as tracer for repeated 
purification of zirconium from accumulated Cb® 
since the minute Cb® has a strong radiation 
intensity. 

The final separation of yttrium from the ex- 
perimental solution of zirconium was made by 
partial and repeated coagulation of zirconium at 
pH, 4.0. 

The common impurities were removed by the 
ammoniacal procedure." 

By direct dilution of the active zirconium 
obtained, as little as 0.06 mg could be used per 
experiment since it had sufficient activity for 
accurate measurement. However, since Zr is 
secured already mixed with stable isotopes the 
ratio of active to stable cannot be increased. 
Therefore, the properties of very dilute solutions 
could be studied only after separation of Zr** 
without carrier, from yttrium bombardment. 

Since the target chemistry in this case is un- 
conventional the properties of zirconium in dilute 
solution will be treated first, returning later to 
the isolation of Zr*®. 


11 Allen, Pool, Kurbatov, and Quill, Phys. Rev. 60, 425 
(1941). 
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Ill. EXPERIMENTAL PROCEDURE 


When ten milliliters of a solution containing 
6X10~- g of zirconium as oxychloride is titrated 
with sodium hydroxide (0.01) opalescence ap- 
pears at a pH as low as 2.8. The solution becomes 
increasingly cloudy as the PH rises then clears at 
a much higher fH than that at which coagulation 
started (Fig. 2). A continuous coagulation process 
is evident during the titration. 

On the basis of these data a pH range well 
below that of quantitative precipitation was 
selected for studies of the behavior of more 
dilute solutions. 

The method used for coagulation and separa- 
tion of precipitate and filtrate was as follows: 
Samples of 20 to 40 ml of known zirconium 
concentration in known acid concentration were 
titrated to a desired pH with standardized 
sodium hydroxide solutions of various strengths 
depending upon the final concentration of zir- 
conium and salt desired in the particular experi- 
ment. These titrations were carried out at 25° 
using a Leeds and Northrup pH meter. Then the 
samples were allowed to stand for an agglomera- 
tion period before removal of the coagulated 
part. 

Since the quality of water is of primary sig- 
nificance in securing reproduceable results the 
following stantlardized conditions were used 
throughout: Portions of the freshly redistilled 
water to be used for a series of experiments were 
electrometrically titrated with sodium hydroxide 
to pH 9 and back titrated with hydrochloric acid 
of known normality to pH 2 to be certain that 
no buffering compounds were present. The initial 
pH of the water was regularly found to_be not 
less than 6.6.* 

The requirement of maintaining PH without 
change during the agglomeration period was met 


* Agglomeration of radio-colloids has been previously 
reported in the literature as attributable entirely to the 
purity of the water used and experiments were described 
using water prepared by centrifuging. In the present work, 
accordingly, centrifuged water was tested in control experi- 
ments with minute concentrations. However, no ob- 
servable effect was found with the centrifuged water as 
compared to the redistilled water maintained free from 
carbon dioxide. (Table II.) In addition, the formation of 
radio-colloids was described previously without reference 
to control of the pH of the solutions or the salts present. 
Consequently the observations made were found to be 
(1931) usive. O. Werner, Zeits. f. physik Chemie A156, 89 
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by performing the experiments in a room free 
from other chemical work using a _ separate 
enclosed space having open vessels of soda-lime 
and provision for air exchange through filters 
with sodium hydroxide..A facilitating circum- 
stance in avoiding carbon dioxide effect is the 
pH range of 3 to 4.5 when the admixture of 
sodium carbonate owing to the initial pH of the 
water constitutes less than 0.1 percent of the 
total salt (NaCl) concentration present in the ex- 
periment. 

With the higher concentrations of zirconium, 
10-5 gA/l and with the higher pH, 4.5, it was 
observed that on standing within the first half 
hour the pH drops to the extent of 0.1 to 0.2. 
On the contrary, no experimentally observable 
change in pH takes place during the agglomera- 
tion period with the very dilute solutions of 
zirconium. The observation of pH lowering on 
standing only with higher concentrations of zir- 
conium was attributed to the fact that the 
agglomeration process requires some time. The 
duration of agglomeration was standardized at 
six hours, this period having been selected upon 
evaluation of the attainment of full agglomera- 
tion after titration of a sample to a definite pH. 
Separations were made after six hours and four- 
teen hours with 2.78X10-* gAZr/l giving the 
following results: 











Percent Percent 
coagulated coagulated 
pH 6 hr. 14 hr. 
+ BS 13.0 12.8 
4.5 26.3 25.4 








Since radioactivity measurements, owing to 
the geometry of samples relative to the G-M 
counter tube or ionization chamber, vary to the 
extent of +3 pércent of the activity measured 
the above figures indicate that extension of 
agglomeration time beyond six hours is not 
required if there is no removal of coagulated 
fraction and no distortion of other conditions. It 
was found that an hour and a half was insufficient 
since the results vary especially at high con- 
centrations. 

After the six-hour agglomeration period the 
solutions were filtered through Schleicher and 
Schull blue ribbon paper, and the deposits on 
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the paper subsequently dissolved by washing 
first with 2N hydrochloric acid and then several 
times with 0.1N hydrochloric acid until the 
papers were inactive. 

The method of removing the coagulated part 
by means of filter paper was chosen after a com- 
parison of data obtained on agglomeration of 
yttrium,* cerium,’ and zirconium.* Filter papers 
were washed before using with salt and hydro- 
chloric acid solutions of »H corresponding to 
that of the experimental series of zirconium. No 
washing, in the analytical sense, can be done on 
the agglomerated fraction on the papers even 
with wash water of corresponding pH and salt 
concentration, therefore, the volume of solution 
left on the filter paper was reduced to an amount 
not exceeding 2 percent by use of 6-cm filter 
papers. 

The effect of filter paper on agglomeration was 
evaluated by arranging consecutive filtrations, 
that is, the filtrate from one funnel immediately 
filling a second paper placed below the first 
funnel. It was found that in minute concentra- 
tions the amount of coagulation on the second 
paper does not exceed 5 percent of the agglomer- 
ated fraction on the first paper. 

While the effect of the filter paper undoubtedly 
cannot be excluded as a part of the process of 
separation, the experimental data suggest that 
the agglomeration of radio-colloids begins with 
adjustment of pH and after a period of time the 
filtration facilitates the removal of coagulated 
particles by localizing their accumulation on 
the paper. 

With concentrations of 10-5 gA/I of zirconium 
and less no coagulated fraction is visible. Conse- 
quently the experimental procedure cannot be 
followed by common analytical methods. The 
complete transfer of agglomerated phase to filter 


(1948) D. Kurbatov and M. L. Pool, Phys. Rev. 65, 61 
44), 

* It was found that agglomeration and settling of radio- 
colloids takes place even in concentrations below 10~-* gA/I 
without filtration. Samples were titrated to the desired pH 
values and allowed to stand for eight days in stoppered 
graduated cylinders. A top portion (about 25 percent) of 
the liquid was slowly pipetted, and the quantity of material 
found in the upper liquid phase was always lower than that 
of the original. However, the concentration of the liquid 
separated by pipette, without filtration, varied con- 
siderably (+16 percent). Similar observations were made 
by re samples after elapse of an agglomeration 
period. 





paper and complete solution from the filter paper 
into a vessel for measuring on the G-M counter 
was controlled by means of activity determina- 
tions. Secondary treatment was made of the 
beakers and filter papers with weighable quanti- 
ties of zirconium. When zero activity was found 
on the beakers and papers the procedure with 
very dilute solutions was considered to have been 
adequately performed. Careful performance of a 
standardized technique which is entirely con- 
trolled by instruments is required instead of the 
visual judgment of the experimenter, when work- 
ing with concentrations below 10-5 gAZr/I. 

The activities of the solutions so obtained and 
the original filtrates were determined either 
(a) directly by means of a Geiger-Mueller 
counter with immersion tube for liquids or 
(b) by means of a Wolf unifilar electrometer. 
In the latter case the solutions were evaporated 
to dryness in Petri dishes of 4-cm diameter and 
the activities measured with an ionization cham- 
ber filled with freon at 20-lb. pressure. 

Both methods have certain disadvantages. 
Zirconium can be kept in solution only in strong 
acid and the very thin glass of the G-M counter 
immersion tube cannot be used long. The elec- 
trometer measurements require preparation of 
solid samples of uniform geometry and there is 
a loss in total activity owing to the depth of the 
measured material. 

In all cases where the ratio of beta- and gamma- 
radiation intensities was required to establish 
the purity of an isotope the electrometer was 
used. 


IV. SOME PROPERTIES OF ZIRCONIUM IN DILUTE 
SOLUTION AT pH 4.0 AND 4.5 


After purification, the oxychloride solutions 
of bombarded zirconium were diluted as pre- 
viously mentioned so that 20 ml of solution had 
an activity of approximately 400 counts per 
minute, measured on the G-M counter.* 

The concentration of such a solution was found 
by usual analysis of undiluted zirconium to be 
1.2310-> gA/l. The more concentrated solu- 
tions were obtained then by adding the desired 
amount of stable zirconium ; the actual quantities 


* The activity of these solutions, expressed in counts, was 
correlated with electrometer measurements expressed in 
units of a constant source of activity. 
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Fic. 3. Coagulation of zirconium at pH, 4.0. 


NaCl 
Curve conc. Final conc. Zr Zr Coagulated 
No. N gA/l gA/l 
la .025 10-8 to 1076 10-7 to 1076 
Ib 025 10-6 to 1074 10-6 to 1075 
Il 005 10-6 to 1074 10° 


(Curve Ib is a continuation of Curve Ia). 


of zirconium used in this series of experiments 
being from 4.110-5 to 1.0X10- g. 

Agglomeration was obtained in the two pH 
regions, 4.0 and 4.5 while the sodium chloride 
concentration was maintained approximately 
0.005N. 

The deviations in measurements of the ac- 
tivities of filtrates and solutions of coagulated 
part did not exceed +2 percent of the total 
counts taken for experiment. 

In Table I the results are given for the co- 
agulate. The log-log plot of coagulated Zr and 
final concentration (Fig. 3, Curve Il) approaches 
a straight line only below 1.2X10~* gAZr/I. If 
the distribution of zirconium between solution 
and coagulate in this region is expressed as 


C/c*=K, 


where c is the amount coagulated per liter of 
solution and C is the final concentration, the 
value of m is 4.1 at 25° with sodium chloride 
concentration ~0.005N and pH 4.0. 

At higher concentrations than 1.210-* gA/I 


AND J. 
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at pH 4.0 and 4.5 the log-log plot of coagulate 
against final concentration does not hold as a 
straight line. However, the log-log plot of co- 
agulate against initial concentration (shown in 
Fig. 4) gives a straight line with exponent 
n=0.48 (pH =4.0) and 0.58 (pH =4.5) and pro- 
portionately constant, k=2.2x10-* and 8.0 
X 10-* respectively. 

Thus the interpretation of zirconium coagula- 
tion as a distribution phenomenon with the 
above dilutions and under the condition that the 
duration of the experiments be maintained con- 
stant is limited to the lower concentrations. The 
time allowed for coagulation was found to be 
significant so, to make results comparable, these 
and following separations of the coagulated 
fraction were performed six hours after titration. 

From the tabulated data it is evident that 
(a) the coagulated fraction increases with dilu- 
tion in this zirconium concentration range, at a 
given salt, pH, and temperature and (b) that 
the fraction of zirconium removed depends on 
PH at this temperature and salt concentration. 

These results were used as a basis for the 
separation of pure Zr** from bombarded yttrium. 


V. TARGET CHEMISTRY OF YTTRIUM AND SEPA- 
RATION OF Zr*® WITHOUT CARRIER 


Yttrium oxide, Hilger’s spectroscopically pure, 
was bombarded with deuterons for four and a 
half hours in the Ohio State University cyclotron. 
The main nuclear reactions are: Y*(d, p)Y%, 
Y89(d, 2n)Zr®®, and Y*°(d, a)Sr*’ and owing to the 


Table I. Coagulation of zirconium from dilute solu- 
tions at 25°. (Final volume approximately 40 ml per 
experiment.) 











Zr coagulated 





Initial conc. NaCl 
Zr gA/l normality gA/l 
104 pH X<10% X104 Percent 
pH 4 
124 3.90 4.9 .103 83 
305 4.01 5.0 162 53 
612 4.02 5.0 .220 36 
1.20 3.97 5.1 252 21 
1.74 4.04 5.2 331 19 
2.78 3.88 5.4 445 16 
pH.4.5 
306 4.43 5.0 211 69 
.606 4.53 5.0 285 47 
1.18 4.44 5.1 413 35 
1.74 4.53 5.2 574 33 
2.76 4.34 i 


.690 25 
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presence of fast neutrons Y*(n, 2n)Y*’, and 
Y*(n, p)Sr®* may occur. 

Since Zr*® disintegrates with emission of posi- 
trons and gamma-radiation and Y°* by emission 
of electrons and no gamma-rays, the extent of 
separation of zirconium from yttrium without the 
addition of stable yttrium can be observed from 
the beta-gamma ratio in Zr*® and the presence 
or absence of gamma-radiation in Y°°.* 

In Fig. 5 is a condensed outline of the pro- 
cedure used in the isolation of Zr**. It may be 
noted that after solution and purification of 
yttrium by the ammonia procedure (not shown 
on chart) yttrium was again purified by solution 
in hydrochloric acid and precipitation with 
sodium hydroxide and finally electrometrically 
titrated to pH 4.0. Because of the high salt con- 
centration and large quantity of yttrium, Zr®® 
was coagulated to the extent of 58.7 percent; 
41.3 percent remaining in solution with the 
yttrium. All the yttrium, within the limits of 
the experimental procedure, was found to be 
left in solution. 

On the right side of the chart are shown the 
results on the purification of zirconium. At the 
lower salt concentration 94.2 percent of zir- 






o pH~ 40 
4 pH~ 45 
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Fic. 4. Coagulation of zirconium in the region of con- 
centrations from 2.8 to 0.12 10~¢ gA/I with NaCl concen- 
tration 0.005N. (Plot against initial concentration of Zr.) 


* Identical geometry of samples and thickness of pre- 
cipitates were easily maintained due to the small quantities 
of material. Radiation contributed by x-rays was included 
in the ratio of beta- and gamma-intensities no separation of 
X-ray intensities being necessary for observation of purity 
of fractions. 


conium was coagulated. The beta-gamma ratio 
remained constant and was not decreased by all 
subsequent coagulations. The purity of the zir- 
conium was also confirmed by following its 
disintegration for several half-lives. 


TARGET CEMESIRY 
Y. + Vv. 
Nuclear adios ‘8. eye mcg on ky Mase 


Active impurities 


Titrated with NoOH 





Solution Y 
VON sat ~~. Precipitate 
+ 41.3% Zr” 58.7 % Zr minute 
Purification of Zr” 
low 
Solution Preaprtate 
zr Zr® 
a 942% 
ACTIVITIES OF FRACTIONS | 
8 radistion 293.0 447 Se 5.16 
y radiction 035 ‘95 O19! 308 
OISTRIBUTION OF Zr & Y in FRACTIONS 
ze 6.3% 35.0% 3.4% $53 
Y 99.59% 0.41% a om 
COMPOSITION OF FRACTIONS 
Zz 2% 73.4% wos 00% 
Yv 99.6% 266% — — 


Fic. 5. Scheme of purification of bombarded yttrium and 
isolation of Zr®*® in pure state. 


During this work an evaluation of the half-life 
of Zr*® was made and found to be slightly longer 
than previously reported owing mainly to better 
removal of yttrium from the zirconium fraction. 
In addition, the nuclear reaction Y**(d, 2n)Zr*® 
was established. 

On the left-hand side of the chart are the 
results on separation of yttrium and zirconium 
with zirconium carrier added in quantity equiva- 
lent to the yttrium. The yttrium left in solution 
in this case was uncontaminated by zirconium 
but the precipitate of zirconium (41.3 percent of 
total) contained 26.6 percent of yttrium. 

At the bottom of the chart are shown relative 
beta- and gamnra-intensities for each fraction and 
the percentages of zirconium and yttrium as 
deduced from these measurements. 

The precision of common gravimetric separa- 
tion of two substances is usually determined by 
a comparison of weights found and known 
weights originally mixed in the sample. In case 
some of each constituent appears mixed with the 
other there is no measure of the purity of the 
fractions. However, with the aid of tracer pro- 
cedure it can be seen that by means of electro- 











214 M. 


TABLE II. Coagulation of zirconium from dilute and minute 
concentrations at 25° and pH =4. 











NaCl Zr coagulated 
Zr conc. normality 
gA/l pH 102 gA/l Percent 

Minute 4.01 2.38 — 69* 
1.541077 4.00 2.41 1.06 1077 69** 
1.54 10-6 4.00 2.46 1.06 10-6 69 
2.69 x 10-6 4.00 2.54 1.67 x 10-6 62 
1.351075 4.00 2.53 6.35 x 10-6 47 
2.69 x 10-5 4.00 ra 1.02 10-5 38 
1.34 10-4 4.01 2.54 3.62 K 10-5 27 
1.54 10-4 4.00 2.40 3.85 X 10-5 25 








* Centrifuged water. ; 
** Water not centrifuged: freshly redistilled. 


metric titration pure zirconium can be secured 
even in minute concentration but part will 
remain in the yttrium fraction. An yttrium frac- 
tion free of zirconium can be obtained if part of 
the yttrium is left in the zirconium. This is 
caused by the overlapping of the pH ranges of 
coagulation of these two elements. 


VI. PROPERTIES OF ZIRCONIUM IN CONCENTRA- 
TIONS BELOW 10-° gA/l 


Using the radioactive Zr**, a series of solutions 
was prepared by adding stable zirconium in 
known quantities, and in order to link these 
experiments with those previously described the 
highest concentration prepared was 1.54X10~4 
gA/l. The lowest concentration of Zr®® was 
approximately 10-" gA/l. This did not contain 
common zirconium and is referred to later as 
minute.* 

The same standardized procedure used on zir- 
conium concentrations above 10-* gA/l was 
followed with the lower concentrations. This in- 
cluded the precautions as to the purity of the 
water and reagents as described above; preven- 
tion of change in pH during the agglomeration 
by maintaining the samples in an enclosed space 
free from carbon dioxide, other gas impurities 
and laboratory dust and especially care with 
regard to the temperature of experiments and 
duration of agglomeration time. The variation 
in temperature allowed in a complete series of 
experiments did not exceed +0.5° from the 25° 
selected. 

* The total yield of purified Zr®® was estimated by com- 


parison of activity with known natural sources to be of the 
order of magnitude of 10-¥ g. 
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Coagulations were carried out at pH 4 with a 
salt concentration of ~0.025N giving the results 
shown in Table II. It can be seen that the frac- 
tion coagulated increases with dilution to a 
certain limit, in this case, 1.54X10-® gA/I. 
Below this concentration the percent of zir- 
conium coagulated remains constant. 

The log-log plot of coagulated and final con- 
centration of zirconium from 1.5X10-* to 2.7 
X10-* gA/I approaches a straight line (Fig. 3, 
Curve Ib). The constant, ”, for this curve at 
PH 4 and 0.025N salt concentrations is 1.5. 

These values can be compared with those 
given previously for pH 4: the fivefold increase 
in salt concentration resulting in less increase in 
coagulation with dilution. 

Within a narrow range of concentrations 2.6 
to 1.54X10-® gA/Il the slope, m, decreases to 
unity and becomes constant and k decreases to 
0.45 (Curve la, Fig. 3). 

It is well known that electrolytes aid in the 
coagulation of colloids and this is true in the 
case of zirconium in concentrations greater than 
10-5 gA/l. In Table III are given the results on 
coagulation at three different pH values and two 
salt concentrations and it can be seen that the 
higher salt results in consistently higher fractions 
coagulated. 

However, with minute concentrations the effect 
of salt is reversed. In Table IV it is shown that 
at a given fH a tenfold increase in salt results 
in a smaller fraction coagulated. 

It should be expected therefore, that at some 
intermediate concentration of zirconium and at 
a certain pH a change in salt concentration would 
be without apparent effect. This is found in the 
zirconium concentration range of 7X10-> and 
10-* gA/l at pH 4 when the salt concentrations 
are as shown: 





Zr conc NaCl Percent NaCl Percent 
gA/l N coag. N coag. 
1.54 X 10-6 .0246 69 503 70 
7.25X 10-5 .025 31 005 31° 








* Data read from crossing of curves Ib and II, Fig. 3. 


The curves in Fig. 6 show that at the same 
minute concentration the fraction coagulated 
within the limits of the experiments is directly 
proportional to the pH, that is, y= (pH—K)/m, 
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r being the coagulated fraction, m, the slope and 
k the intercept. The constants are m=3.15 for 
both curves and k=1.85 and 1.58 for 0.05 and 
0.005N salt respectively. 


The effect of ~H on coagulation of zirconium. 


throughout the range of concentrations studied, 
has been observed to be in the same direction in 
all cases, that is, increase in pH increase coagula- 
tion, salt being constant. 

The basic observation on the coagulation of 
zirconium from very dilute solutions which is in 
line with the results previously observed with 
yttrium, is that this hydrous oxide does not 
exhibit a constant solubility value. 

This has been confirmed by experiments ex- 
tended for eight days. Two different concentra- 
tions of zirconium were coagulated at three 
different pH values. The results are shown in the 
left column of Table V. After eight days the pH 
was redetermined and a second coagulated part 


TABLE III. Effect of salt concentration and pH on 
dilute solutions of zirconium. (Final volume ~45 ml/ 
experiment.) 











Zrconc. NaCl Zr conc. 
gA/l conc. Percent gA/I NaCl Percent 
x104 N X108 pu coag. <x104 NxX<102 pH coag. 
2.86 3.43 13 1.91 3.47 17 


5.3 2.55 
246 SA 33 % 2.50 1.88 4.00 28 
5.5 447 25 2.44 182 443 98 








TABLE IV. Effect of salt and pH on coagulation of minute 
concentrations of zirconium. 





= 








NaCl conc. Percent NaCl conc. Percent 
N X108 pH coag. N X10? pu coag. 
4.57 2.50 29 5.01 2.50 20 
5.53 3.01 47 5.35 3.00 35 
5.85 3.51 60 5.40 3.52 53 
5.93 4.00 77 5.43 4.02 63 

5.43 4.25 76 








TABLE V. Two successive coagulations of zirconium. 








Second coagulation from 


First coagulation 
filtrate. Separation after 8 days 


separation after 6 hours 





Conc. Zr Per- NaCl Conc. Zr Per- NaCl 
gA/l cent conc. gA/l cent conc. 
104 pH coag. N X103 x<104 pH coag. N X10? 

622 3.48 24 4.9 1.29 sae kt i. 

612 4.02 36 5.0 .950 342 31 186 

606 4.53 47 5.0 Pi (nt aS Me eg 
1.80 352 >». oa 4.07 3.34 .. :o 
1.74 4.04 19 5.2 3.46 3.73 14 1.1 
1.74 53° 23° G2 2.86(2) wee FF) Be 
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Fic. 6. Effect of salt and pH on a minute concentration of 
zirconium at 25°. 


removed from which it can be seen that at essen- 
tially identical salt concentration and pH (pH 
— 3.97: NaCl—.011N) a three and a half-fold in- 
crease in concentration of zirconium in solution 
has negligible effect on coagulation.* 

The experiments described may be partially 
interpreted if the phenomenon of distribution of 
zirconium between solution and coagulate be 
accepted as a Nernst distribution. The law holds 
well for concentrations below 10-* gA/I and up 
to 1X10~* gA/I if the salt concentration is as 
high as 2.5X10-2N. 

It has been pointed out that in heterogeneous 
systems the distribution. law holds well only for 
dilute conditions. With the lowest concentrations 
of zirconium studied, it was found that the 
amount coagulated is directly proportional to the 
amount left in solution. In the light of classical 
concepts of distribution this would mean only 
one molecular species of zirconium in such dilute 
solutions under the conditions specified. 

In the region of 10-* to 10-* gAZr/I in 
which the log-log plot of zirconium in solution 
against zirconium coagulated closely approaches 
a straight line the exponent 1, is 1.5 with 0.025N 
salt and 4.1 with 0.005N salt. This may indicate 
that the zirconium left in solution exists in more 
than one form in this range or that the coagu- 
lated part is more associated. 

The log-log plot of coagulated Zr against final 
concentration of Zr does not give a straight line 
for concentrations >10-* gA/I. This indicates a 
more rapid change in the composition of zir- 


*If the solubility were constant and these solutions 
saturated the coagulated fraction, for example, from the 
solution numbered (2) should be 92 percent instead of 77 
percent if the solution numbered (1) were considered 
saturated. 
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conium in solution or coagulate than that ob- 
served in the range immediately above 10-6 gA/I. 

It has been established that coagulation of Zr 
will take place at any concentration under the 
conditions described. Accordingly, it was shown 
above that for consecutive coagulations a second 
addition of hydroxyl ions is not required. The 
process proceeds at somewhat lower pH values 
but extremely slowly. 

The observation that the fraction coagulated 
increases with dilution to some limit is well 
supported by the effect of hydroxyl ions on 
hydrolysis. 

Below 10-* gA/I the coagulation of zirconium 
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reaches a maximum and the fraction remains 
constant consequently the separation of Zr®® in 
pure state, in quantities of the order of 10-'° g, 
was possible. 
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The low velocity scattering of H+, H2+, and H;* in methane was measured at voltages be- 
tween about 1 and 150. Neutralization is found in all three cases, being much greater for H2* 
and H* than for H;*. The elastic scattering is also relatively high compared to hydrogen and 
helium. It was found that there was an effect of pressure at constant voltage on the experimental 
determination of effective cross section for elastic scattering. This has been explained on the 
basis of multiple scattering. Approximate potential laws for the interactions have been evaluated. 


HE experience gained in previous experi- 
mental! work on the measurement of low 
velocity ion scattering and the techniques de- 
veloped in this work enabled-the measurements 
of the scattering of low velocity hydrogen ions in 
hydrocarbons to be approached with confidence. 
In particular, the techniques of the determination 
of neutralization’ and of the removal of 
carbonaceous films! provided procedures which 
enabled such measurements to be made. Methane 
was the choice for the first hydrocarbon to be 
1 (a) A. S. Russell, C. M. Fontana, and J. H. Simons, J. 
Chem. Phys. 9, 381 (1941). (b) J. H. Simons, H. T. Francis, 
C. M. Fontana, and S. R. Jackson, Rev. Sci. Inst. 13, 419 
(1942). (c) J. H. Simons, C. M. Fontana, E. E. Muschlitz, 
Jr., and S. R. Jackson, J. Chem. Phys. 11, 307 (1943). 
(d) J. H. Simons, C. M. Fontana, H. T. Francis, and L. G. 
Unger, J. Chem. Phys. 11, 312 (1943). (e) J. H. Simons, 
H. T. Francis, E. E. Muschlitz, Jr., and G. C. Fryburg, J. 
Chem. Phys. 11, 316 (1943). (f) J. H. Simons, E. E. 


Muschlitz, Jr., and L. G. Unger, J. Chem. Phys. 11, 322 
(1943). 


studied because of its simplicity and gaseous 
nature. 

In the same apparatus and using the same ex- 
perimental methods and analytical correlations 
as have been previously described,' the elastic 
scattering and neutralization in methane of all 
three hydrogen ions, H+, Het, and H3;+t, have 
been measured from approximately 1 to 150 
volts. The formation of carbonaceous films from 
the methane required their frequent removal 
(between every several experiments) by the 
method of a discharge in oxygen. The results are 
given in Table I and shown graphically in 
Figs. 1-3. 

In the text, Rr represents the ratio of the 
unscattered current to the total ion current 
entering the scattering chamber; and Rs repre- 
sents the ratio of the unscattered current to the 
total current minus that deflected by neutraliza- 
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SCATTERING OF HYDROGEN 
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Fic. 1. Scattering of H* in methane. 











TABLE I. Scattering of hydrogen ions in methane. 
Ht H2* H;* 
Ve- Ve- Ve- 
locity locity locity 
in ; in in 
Volts aT a,” Volts aT a, Volts aT as 
156.0 113.0 43 125.0 110.8 42 152.0 35.0 24 
130.0 117.0 47 125.0 109.6 41 130.0 35.1 25 
128.0 116.0 40.5 124.0 112.8 44 128.0 36.4 24.4 
110.0 120.0 47 115.0 117.0 46 125.0 34.0 25 
102.0 123.0 44 99.0 120.0 47 116.0 38.2 25.2 
96.0 126.0 51 96.0 118.8 46 108.0 36.6 28 
77.5 132.0 55 92.0 121.0 47 103.0 39.5 26.7 
72.0 135.5 51 89.0 123.0 48 101.0 39.2 27.7 
58.0 144.0 60 78.0 128.0 50 80.3 41.2 27 
52.5 151.0 62 78.0 125.1 49 78.0 40.9 31.5 
52.0 153.0 55 71.0 133.0 52 74.0 42.6 30.1 
42.8 161.0 59 62.0 134.0 53 66.0 43.5 32.3 
40.9 160.0 65 S75 1366 33 63.5 46.0 31 
31.0 173.0 69 50.55 141.0 56 56.0 45.0 34.5 
24.6 180.0 70 45.0 145.5 58 49.0 49.7 35.5 
23.0 187.0 70 375 1519 62 49.0 50.4 33.3 
19.5 188.0 75 28.5 162.0 64 47.0 52.6 37 
17.0 193.0 73 25.0 165.0 66 44.0 50.8 35 
15.0 198.0 77 25.0 169.0 62 34.5 55.1 38 
11.6 214.0 78 18.0 181.0 67 32.0 57.2 42 
8.7 233.0 86 11.5 188.0 70 30.5 58.8 43 
7.8 220.0 91 8.1 209.5 83 28.0 62.9 42 
44 233.0 99 5.1 220.0 82 26.5 62.1 42.4 
1.1 289.0 128 3.0 248.0 96 20.5 68.2 48 
13.6 88.3 61.5 
13.6 87.2 61.5 
12.0 90.9 62 
11.5 86.8 60 
11.5 95.3 65 
9.5 98.7 69 
8.5 108.3 78 
7.8 113.0 82 
3.9 157.0 112 
3.0 160.0 111 
1.5 185.0 148 











Gas pressures of the order of magnitude of 1 X10-3 mm were used for 
H* and Het and 2 to 3 X10~ for Hs+. Empty tube pressure was 1 X1075 
mm or less. Temperature was 25 to 30°C. ay® =a7 —as°. 
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Fic. 2. Scattering of H2* in methane. 





200 -— — 














0 $0 Win Volls 100 10 


Fic. 3. Scattering of H;* in methane. 
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Fic. 4. Effect of pressure on elastic scattering. 


For scattering in methane of 56 volt H+(©), 56 volt H2t(O), 125 volt 
H2*(Q), 56 volt Hs*(Q ), and 125 volt Hs*(@). 


tion. Rr® and Rs° are the same ratios taken at 
some very low pressure in the measuring tube to 
correct for tube alinement, etc. ay, as, and a7 are 
the experimental effective cross sections for 
neutralization, elastic scattering, and total scat- 
tering, respectively, obtained without correction 
from the measured ratios. as® and ay® are the 
corrected values of as and ay. ay°=a7—as°. 

A further difficulty was encountered in these 
measurements. This had been theoretically an- 
ticipated previously but not experimentally ex- 
perienced. There is a possibility of an ion that has 
been deflected from the primary beam by an 
elastic scattering collision to suffer a further 
collision before reaching the wall of the scattering 
chamber. If this second collision is an elastic one, 
the effect on the measurement will be nil, as the 
chance of deflecting this scattered ion into the 
collecting chamber is very small. The chance of 
the elastic collision changing its direction suff- 
‘ ciently to direct it to the scattering chamber lid 
is also extremely small, and these are the only 
two possibilities that would influence the meas- 
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urement. However, if the second collision is an 
inelastic one (neutralization), then the ion already 
scattered: will be collected on the scattering 
chamber lid and so instead of being counted as 
elastically scattered, it is counted as neutralized. 
The total scattering would be almost unaffected 
by this process but the relative amounts of elastic 
scattering and neutralization would be affected. 
This effect would be most marked for substances 
having relatively large values of neutralization, 
and it would also be a function of the pressure, 
increasing with increasing pressure. As the pres- 
sure increased, the experimentally determined 
elastic scattering would be too low and the 
neutralization too high. To obtain more accurate 
measurements of elastic scattering for any sub- 
stance experiencing considerable neutralization, 
measurements at any one voltage could be made 
over a range of pressures and these extrapolated 
to zero pressure. As 


(Rs°/Rs)—1 
as= aT; 
(Rr°/Rr)—-1 





TABLE II. Variation of elastic scattering with pressure. 











P X105 Fract. 
Particle (mm) Scatt. aT as as” ay? 
125 volt 32.53 0.645 3o4 226 24 11.1 
29.53 0.655 Jon 22.4 23.9 10.3 
H;* 28.43 0.658 | ee) ee 10 
21.20 0.716 34.0 244 25 9 
15.38 0.722 33.0 23.8 24.5 8.5 
Extrap 0.795 33.9 27 6.9 
56 volt 57.83 0.592 44.9 266 31 13.9 
34.33 0.710 44.1 31.3 53.2 10.4 
H;* 20.29 = 0.717 9 #5 33 10.9 
11.95 0.764 43.9 33.5 34.2 9.7 
6.47 0.785 43.2 33.9 34.2 9.0 
Extrap 0.790 44.0 34.7 9.3 
125 volt 28.72 0.269 1088 29.3 42.5 66.3 
15.38 0.312 1096 34.2 42 67.6 
H.* 14.83 0.314 110.8 34.7 42.6 68.2 
8.86 0.330 109.6 36.2 41 68.6 
Extrap 0.360 109.7 39.5 70.2 
56 volt 22.89 0.282 1500 42.3 63 87 
15.66 0.300 150.3 45.1 60 90' 
H.* 11.24 0315 1520 47.9 59 93 
7.4 0.342 1530 52.33 60 93 
4.39 0.344 1540 53.0 58 96 
Extrap 0.365 154.0 56 98 
S6vokt 2343 O241 15759 318 3B 98 
17.69 0.263 1560 41.0 58 98 
Ht 7.8 0.320 161.0 51.5 60 101 
3.644 0.328 160.0 52.5 57 103 
Extrap 0.360 160.0 57.7 102.3 
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Fic. 5. H* and H;* in methane. 


and 


(Rs°/Rs)—1 
(Rr°/Rr)—1 





is obviously the fraction of total scattering that 
is scattered elastically, a plot of this fraction 
against pressure at constant voltage should be 
extrapolatable to zero pressure and this extrapo- 
lated value used to obtain as® from ar. Such a 
plot is shown in Fig. 4. In Table II the experi- 
mental values are given, and values of a7 and as, 
at constant voltage but at varying pressures are 
given. It is seen that a7 is practically independent 
of pressure. as decreases with pressure. 

In collecting the experimental data for this 
study the large contribution of this pressure effect 
was not appreciated until much of the data had 
been collected and the pressure effect not experi- 
mentally determined at every voltage. In the 
following paper on the scattering of hydrogen ions 


in ethylene and propylene an effect of pressure on 


the total cross-sectional area is described. This 
was not experienced in this work, probably be- 
cause of somewhat different experimental condi- 
tions and because the total amount of scattering 
was smaller. An analytical treatment is given in 
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this following paper for the calculation of the 
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Fic. 6. H»* scattered in methane. 


cross section for elastic scattering corrected for 
the effect of pressure, as°, which does not require 
measurements at a variety of pressures at the 
same velocity. We have used this method of 
calculating as®. In Table I] we compare these 
values with the value of as° obtained by extrapo- 
lating the experimental values to zero pressure. 
A satisfactory agreement is found. 

The determination of a potential law from the 
elastic scattering curves for these experiments is 
less sure than for previous work for three reasons. 
Because of the greater total scattering power of 
methane, the pressure used in the scattering 
chamber is less and is thus determined with less 
precision. The electrical measurements are made 
with greater difficulty caused by the greater 
tendency to form carbonaceous films. The error 
arising from these two causes is of the order of 
4 percent. In addition, the effect of pressure at 
constant voltage on the elastic scattering intro- 
duces some uncertainty. 

In order to calculate potential laws for the 
elastic scattering, log as° is plotted against log W, 
where W is the velocity in volts. Figure 5 shows 
such plots for H+ and H;* and Fig. 6 for H2*. In 
these figures the straight lines are obtained by 
using values of the coordinates of the smooth 
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curves of as® in Figs. 1-3. The points are calcu- 

lated from the experimentak values. Potential 

laws calculated from these lines are as follows. 
For H+ in methane 


—7.0X10-® with K in volt cm® and 7 in cm 
r? between 5 and 40 volts (2.96A 
and 2.38A) ; 


—7.8X10-*! with K in volt cm’? and ¢ in cm 




















wie ad between 40 and 150 volts (2.38A 
and 1.97A). 
For H2* in methane 
—1.2X10-* with K in volt cm" and 7 in cm 
- yi0 between 5 and 40 volts (2.8A 
and 2.29A) ; 
ae 1.5X10-® with K in volt cm® and 7 in cm 
? ~ r between 40 and 130 volts (2.29A 
and 1.97A). 
For H;* in methane 
ad, —2.8X10-*! with K in volt cm‘ and 7 in cm 
yf between 5 and 22 volts (2.96A 
and 2.07A) 
Be | —7.8X10-** with K in volt cm® and ¢ in cm 


r® between 22 and 120 volts 


(2.07A and 1.5A). 


The constants are not corrected for mass ratio as 
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shown in previous work! with the exception of 
the one which is starred. 

This starred potential function may be com- 
pared to the potential function which would exist 
if a simple polarization force acted between the 
H;+ ion and the methane molecule. R. Sanger? 
has given a value for the dielectric constant of 
methane expressed as (e—1), where e€ is the 
dielectric constant, of 84.2 10-5 at 291.7A. From 
this, the polarizability may be calculated as- 
suming only that methane obeys the ideal gas 
law at this temperature and at one atmosphere 
pressure. The potential function obtained from 
the polarizability is , 


— 1.90 10-*! 





with K in volt cm‘ and? in cm. 
4 
r 


This is in reasonable agreement with the experi- 
mental potential function given above expressing 
the interaction between H;*+ and the methane 
molecule at large distances. 

There are several interesting facts uncovered 
by these experiments: the exhibition of the 
neutralization phenomenon by H;+t—this is the 
first case in which H;+ has undergone neutraliza- 
tion; and the large percentage of neutralization 
displayed by H+ and H2+—about 70 percent of 
total scattering. This latter indicates a strong 
tendency on the part of methane readily to lose 
an electron forming the positively charged aggre- 
gate CH,* in-preference to attracting a positive 
ion to form the aggregates CH;*+ and CH,¢t. The 
high values for the exponent in the potential laws 
are also of interest. 


2 R. Sanger, Physik. Zeits. 27, 556 (1926). 
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The low velocity scattering of H*+, H2*+, and H;* in both ethylene and propylene has been 
measured at ion velocities of less than 150 volts. Both elastic scattering and neutralization have 
been obtained over a velocity range. An effect of gas pressure at constant voltage on the total 
scattering was experienced. This is briefly discussed as is also the change of the relative amounts 
of elastic scattering and neutralization with pressure at constant voltage. The amount of 
neutralization compared to elastic scattering is large in all six cases. The elastic scattering in 
ethylene is somewhat similar to that in methane whereas that in propylene is somewhat greater. 





. a previous paper! experiments on the elastic 

scattering and neutralization of hydrogen ions 
in methane were reported. In exploring the possi- 
ble use of this technique for organic chemical 
substances, the olefins ethylene and propylene 
were selected for the present study. The previ- 
ously reported techniques were employed, and 
it was found possible to make satisfactory 
measurements. 


EFFECT OF PRESSURE ON THE TOTAL CROSS- 
SECTIONAL AREA 


Inability to check the total cross-sectional 
areas in experiments at different total gas pres- 
sures led to an investigation of the total cross- 
sectional area as a function of pressure. It was 
found that a plot of ar (total cross-sectional area) 
against pressure at constant velocity gave ap- 
proximately a straight line. At a lower velocity 
where ar is larger, the variation of ar with pres- 
sure is much greater. A quadratic relationship 
was found to hold between the slopes of the above- 
mentioned lines and the value of a7 extrapolated 
to zero pressure, a7’. 


—dar/dp=K(ar)’. (1) 


Integrated between pressure limits, » and zero, 
this becomes 


ap’=ar+K(ar®)p. (2) 


For convenience and without serious loss, (a7r°)? 
is replaced by ar” to give 


ar’=ar(1+Karp). (3) 


1]. H. Simons, and G. C. Fryburg, J. Chem. Phys. 11, 
132 (1943). 
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The value of K in Eq. (3) appeared to be ap- 
proximately the same for the three different ions 
in both ethylene and propylene and also for some 
of the data previously reported on the scattering 
in water vapor,’ when reexamined in the light of 
this effect. Perhaps, therefore, K is a property of 
the instrument and a function of the constants of 
the apparatus. The value of K thus found is 0.35. 

In all previously reported work on low velocity 
ion scattering, except that of water, this correc- 
tion is smaller than the experimental error. Its 
value depends upon the product arp and where 
this is small, the total cross section is inde- 
pendent of the pressure within the precision of 
the measurements. In the scattering in water 
vapor the value of arp is larger, but the total 
precision of the measurements is much less due to 
the difficulty of pressure measurements. This 
correction applied to the scattering in water 
vapor is probably not greater than the experi- 
mental error. 


EFFECT OF PRESSURE ON THE FRACTIONAL 
AMOUNT OF THE TOTAL CROSS-SECTIONAL 
AREA TO BE CREDITED TO 
ELASTIC SCATTERING 


A previous paper! described an effect of pres- 
sure at constant voltage on the effective cross- 
sectional area for elastic scattering. An ion in the 
beam in the scattering chamber is counted as not 
suffering a collision of any kind if it is caught in 
the collecting chamber. If it reaches the walls of 
the scattering chamber, it is initially recorded as 


2]. H. Simons, H. T. Francis, E. E. Muschlitz, Jr., and 
G. C. Fryburg, J. Chem. Phys. 11, 316 (1943). 
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suffering an elastic collision; and if its resultant 
of an electron exchange reaches the lid of this 
chamber under the influence of the retarding 
field, it is initially counted as having been 
neutralized. After the primary act of neutraliza- 
tion, a further electron exchange or elastic colli- 
sion will not effect the determination. Also, after 
a primary act of elastic scattering, a further 
elastic scattering will effect the measurements 
only under the condition that this second collision 
diverts the ion either into the collecting chamber 
or to the lid. However, if after an initial elastic 
scattering, the ion suffers neutralization, the 
measurements will need correction as it will now 
be counted as neutralized rather than elastically 
scattered. For experiments where neutralization 
is relatively large this correction becomes 
appreciable. 

If Iv, is the amount of ion current which is 
scattered and subsequently neutralized, J, is the 
total current scattered, J,’ is the uncorrected 
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Fic. 1. Scattering of H* in ethylene. 
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scattered current, x is the distance the scattered 
ion travels before neutralization, d is the average 
.distance traveled by a scattered ion before 
reaching the scattering cylinder, and p is the 
pressure, then 


dIns=TIan%pdx, (4) 
I,’ =I,—Ins, (5) 

by integration between 0 and d. 
In I,/(I,—Ins) = pawd. (6) 


It can be seen that the ratio of the cross-sectional 
area for elastic to total scattering is now modified 
to 


a,®/ar®=[(1/R,—1)/(1/Rr—1) ]J./7’, (7) 
or 
log a,°/ar®=log [(1/R.—1)/(1/Rr—1) ] 
+ (d/2.302)an°p. (8) 
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Fic. 2. Scattering of H2* in ethylene. 








SCATTERING OF HYDROGEN 


IONS IN PROPYLENE 


TABLE I. Scattering of hydrogen ions in ethylene. 
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TABLE II. Scattering of hydrogen ions in propylene. 
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A corrected value of a, can be obtained from 
Eq. (8) if log [(1/R.—1)/(1/Rr—1) ] is known at 
zero pressure. An extrapolated linear plot be- 
tween these variables will give this for measure- 
ments at the same velocity but at different 
pressures. A value of day is obtained from the 


slope of this line. However, if d is known, a,° may 
be obtained from a single measurement from 
Eq. (8) by means of successive approximations. 
From a large number of measurements an aver- 
age value of d was found to be 2.0 cm. As this is 
about one-half the length of the scattering 
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Fic. 3. Scattering of H;* in ethylene. 


chamber, it is a reasonable value. The reported 
values of a, are obtained using this value of d 
and Eq. (8). 


MATERIALS AND PROCEDURES 


Both the ethylene and the propylene used were 
obtained from commercial cylinders and purified 
by repeated simple distillations. Check determi- 
nations were made with different samples. The 
experimental procedures used were the same as 
has been previously described. 


RESULTS 


The results of the experiments on the scattering 
of H+, H.+, and H;+ in ethylene are given in 
Table I, and shown in Figs. 1, 2, and 3; and those 
for scattering in propylene are given in Table I] 
and shown in Figs. 4, 5, and 6. Electrical measure- 
ments were made down to 0.5 volt, but below 2 
volts the precision of measurement is question- 
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Fic. 4. Scattering of H* in propylene. 


able due to the pressure effects and errors in the 
pressure measurement. Gas pressures of the order 
of magnitude of 1X10-* mm were used. Empty 
tube pressure was 2X10-5 or less. Temperature 
was 25 to 30°C. ay®=ar°—a,°. 


DISCUSSION 


It was more difficult to obtain results of pre- 
cision comparable to those obtained with hydro- 
gen and helium. Although the use of hydrocarbon 
gas increased the experimental difficulties due to 
the formation of carbonaceous films, the chief 
reason for the lowered precision results from the 
fact that the total scattering is very much higher; 
and to make measurements in the same measuring 
tube, a lower pressure had to be employed. The 
poorer consistency of the measurements resulted 
from the difficulty of measuring this lower 
pressure. 

The nature and the magnitude of the inter- 
action of hydrogen ions with these olefins is a 
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matter of considerable inicrest. It is to be 
observed that a relatively large amount of 
neutralization was obtained. The elastic scatter- 
ing is also found to be of considerable magnitude, 
that for ethylene being similar to that found 
previously for methane and that for propylene 
being somewhat higher. Despite the lower pre- 
cision of the elastic scattering data, it was 
thought desirable to explore the force laws that 
could be evaluated. In Figs. 7, 8, and 9 are given 
the logarithmic plots for’ the evaluation of the 
force laws between ethylene and Ht, Hg", 
and H;+, respectively. The full lines shown in 
these plots are the lines used for the potential 
functions for these interactions, as shown in 
Table III. 

A potential function based upon the interaction 
of a charge with a polarizable molecule contains 
the inverse fourth power of the distance. From 
the polarizability of the molecule, the constant 
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IONS IN PROPYLENE 225 
and the function can be obtained. Debye® has 
given the polarization of ethylene as 9.4 cc per 
mol. The potential function based on this simple 
interaction with ethylene is 


V=—26.8/r', 


with K in volt (Angstroms)! and r in Angstroms. 
It is seen, when the straight lines that will give 
an inverse fourth-power law are drawn on Figs. 7, 
8, and 9 and the constants corrected for mass 
ratios, that these constants found for Ht, Het, 
and H3¢+ are 31.3, 21.8, and 56, respectively. The 
first two are in reasonable agreement with the 
theoretical polarizability law. 

McAlpine and Smyth‘ have given the dipole 
moment of propylene to be 0.35 D and the 
polarization to be 15.7 cc per mole. The simple 
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Fic. 6. Scattering of H;* in propylene. 


3’ Peter Debye, Polar Molecules (Chemical Publishing 
Company, Inc., New York, 1929). 

4K. B. McAlpine and C. P. Smyth, J. Am, Chem, Soc., 
55, 453 (1933). 
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Fic. 7. H* in ethylene. 


potential law based upon the interaction of a 
charge on a molecule of this polarizability and 
dipole moment is 


V = —44.8/r4*—1.05/r?. 


The points in Figs. 10, 11, and 12 cannot be 
forced into a potential law of low exponent. 
However, a line has been drawn on these graphs 
which fits the above law and which crosses the 
path of the points. It is seen that the experimental 
evidence is hardly consistent with this simple 
theory. 

In Figs. 7 and 9 the experimental points appear 
to follow an ‘“S” type of curve. This is an 
interesting observation relative to the variation 
of the force law with distance between particles. 
If the distance decreases from relatively large 
values, the exponent in the attractive law of force 
would be expected to decrease. However, as the 
distance becomes smaller, repulsive forces will 
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TABLE III. Summary of potential functions. 








Range Potential function 





Ethylene 
Theoretical interaction 
of charge and polar- 
izable molecule. 7= —26.8/r' 
3-19 v (3.36—2.45A) = —241/r6 
25-100 v (2.18—2.00A) V=—4.8xX 108/r* 
6-25 v (3.08—2.18A) V=—38.2/ri=[—31.3/r']}* 


Below 5 v (>3.00A) V = —26.8/r!=[—21.8/r*}* 
3-100 v (3.18—2.52A) V=—9.6X107/r'® 


3-25 v (4.68 —2.30A) = —31.3/r' 
25-150 v (2.30—1.87A) V=—990/r8 
3-80 v (4.68—1.90A) V=-—69/r! 


Propylene 

Theoretical interaction 

of charge with polar- 

izable molecule con- 

taining a dipole. 

H+ 6-120 v (2.92—2.12A) 
H.* 5-110 v (3.1—2.8A) 

H;+ 4-100 v (3.66—2.16A) 


V = —44.8/r'—1.05/r 
V=—105/r2 
V=—7.4X 10/7 
V=—537/r6 








vy is given in angstrom units, and all values of K are uncorrected for 
mass ratio except those marked *. These are corrected by the method 
described in J. H. Simons, E. E. Muschlitz and L. G. Unger, J. Chem, 
Phys. 11, 322 (1943). 
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come into play and the laws expressing these have 
large values of the exponent. In these scattering 
measurements, attraction and repulsion are 
indistinguishable so that plots such as shown in 
Figs. 7 and 9 might be expected to form an ‘S”’ 
curve. These are the first cases showing an in- 
creasing value of the exponent ‘‘n”’ with decreasing 
distance. In most previous experiments the value 
of ‘‘n’”’ decreased with decreasing distances, but 
these experiments were for molecules of much 
smaller size. If attempts were made to evaluate 
the energy of combination of, for example, H+ 
and ethylene, it seems that the potential function 
with high values of the exponent at the smaller 
distances should not be employed, for if these are 
due to repulsion, they are recorded at greater 
distances than the equilibrium distance. These 
appear in the experimental results only from the 
statistical effect of the molecules being oriented 
at all angles in space. If, for example, in ethylene 
the equilibrium position of an added proton is 
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IONS IN PROPYLENE 227 
between the two carbon atoms, the average effect 
of the repulsion of the ends of the molecule on the 
ion would be felt by an approaching proton at 
greater distances than the equilibrium distance. 
The reverse calculation of the above has been 
made in this case. Baughan, Evans, and Polanyi® 
have calculated the proton affinity of ethylene to 
be 174 Kcal. Using the sixth-power potential 
function, an equilibrium distance of 1.50A is 
obtained. This compares with 1.34A which is the - 
distance the nucleus of the hydrogen atom of 
ethane is from the center of gravity of the mole- 
cule. It could be reasonably expected that the 
distances in the C,H,* ion might be greater than 
in the neutral molecular ethane. 

If we picture an ion passing an unsymmetrical 
rotating molecule, we can anticipate two alterna- 
tive situations to develop. If the interaction is 
dependent chiefly upon the electrostatic attrac- 
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5 E. C. Baughan, M. P. Evans, and M. Polanyi, Trans. 
Faraday Soc. 37, 377 (1941). 
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tion between the charge of ‘the ion and the 
induced dipole of the molecule and if the 
polarizability of the molecule is isotropic, the 
rotations should be unaffected and the repulsive 
forces experienced upon close approach should be 
felt on nearing the greatest radius of the molecule. 

If, however, forces. of attraction for the ion 
already present in the molecule or induced by its 
approach are not spherically symmetrical, such 
as would be expected for the chemical binding of 
the two particles, a different situation would 
result. As the ion approached the molecule more 
‘closely, the rotation would be converted into an 
oscillation of the same energy and that direction 
of the molecule having the greatest attraction for 
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the ion would be uniformly presented to it. The 
distance of closest approach would not be a 
statistical effect but would be that distance given 
by the properties of the molecule along this 
direction. 

If this argument is applied to the interaction of 
H+ with C;H,, it would be expected that a small 
distance of close approach would be obtained 
before the potential law had a high value of the 
exponent. It might also be expected that the 
value of ‘‘n’’ would decrease to below 2 as was 
found in the interaction of H+ with both H». 
and He. The results obtained indicate that this is 
not the nature of the interaction for H+ with 
ethylene and propylene. 
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The general theory of the adsorption of solutes from solutions which are treated by a process 
similar to that employed in chromatography is presented. By assuming a chemical equilibrium 
between adsorbed material and the solution, the adsorption equations can be solved explicitly. 
A detailed discussion of the two-component problem is given; the many-component problem 
and the extension of the theory to more general adsorption isotherms are considered briefly. 





INTRODUCTION 


HE problem of the adsorption of solutes 

from a solution which is passed through a 
tube containing an adsorbing material has been 
considered from the theoretical viewpoint by 
several authors.'~* The equations which apply to 
this problem are the set of simultaneous partial 
differential equations 


(dc;/dl) + (0g;/dv)=0, i=1,2,---m, -(1) 


where c; is the concentration of the 7’th solute in 
the solution at a distance / along the tube after a 
volume v has passed this point, q; is the concen- 
tration of the corresponding component on the 
adsorbing material, and m is the number of 
solutes which are being adsorbed. Wilson! showed 
that in certain instances g; may be a discontinu- 
ous function of /. DeVault? solved the problem 
for one solute completely, in terms of a general 
adsorption isotherm, and discussed the problem 
of multiple adsorption qualitatively. DeVault’s 
investigation of the one-component problem 
showed that gq will be a discontinuous function of 
l only if the solute is adsorbed strongly. Weiss* 
solved the one-component problem for several 
different types of adsorption isotherms, with 
results similar to those of DeVault. 

In this article the general solution of the 
problem of multiple adsorption is given, as- 
suming that the interaction between dissolved 
and adsorbed material is governed by the ordi- 
nary relations of chemical equilibrium. This as- 
sumption is probably valid for ionic adsorption; 
for example, in the replacement of the hydrogen 

1 J. N. Wilson, J. Am. Chem. Soc. 62, 1583 (1940). 


2D. DeVault, J. Am. Chem. Soc. 65, 532 (1943). 
3]. Weiss, J. Chem. Soc. London 297 (1943). 


ion of hydrogen zeolite by an alkali ion.‘ In addi- 
tion, the solution of this particular problem will 
indicate the results to be expected if the adsorp- 
tion is governed by more general adsorption 
isotherms, since the results depend primarily on 
the curvature of the isotherm and not on its exact 
analytical form. The effects of diffusion are not 
considered. 


ONE COMPONENT 


The equations in this section could be obtained 
from the general equation given by DeVault.’ 
Because we are dealing here with a particular 
type of adsorption it will be convenient to state 
the results explicitly. The solution of the one- 
component problem is valuable as a guide to the 
more general solution. We will visualize the 
process as the replacement of the hydrogen ion of 
hydrogen zeolite by an alkali ion, although of 
course the results apply to any adsorption process 
which can be regarded as a chemical equilibrium. 
The reaction involved could be 


HZ +MX=MZ+ HX (2) 
(P—q) (c) (q) (ce) 


where P is the total replaceable hydrogen per unit 
length of tube, c is the concentration of dissolved 
M, q is the concentration of adsorbed M, and c°® 
is the original concentration of dissolved M. q and 
c will be functions of both / and v. We assume 
that the relation between g and c is given by the 
equilibrium equation 


K=4q(¢°—c)/c(P—q). (3) 


4 Private communication from Dr. F. R. Duke, Frick 
Chemical Laboratory, Princeton University. An experi- 
mental article on this subject will be published by Dr. 
Duke in the J. Am. Chem. Soc. 
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Fic. 1. Distribution of adsorbed material: one component. 


Using (3) to eliminate c from Eq. (1), we obtain 
oq 0q 
KPo—+[KP+(1-—K)q}*—=0, (4) 
al Ov 


as the differential equation which determines q 
as a function of / and v. This differential equation 
can be solved by the method of characteristics.® 
The general solution is 


q=~{[KP+(1—K)q}1—KPc'v}, (5) 


where ¢ is any arbitrary function of the argument 
in brackets. This general solution may be verified 
by direct substitution in Eq. (4). The boundary 
conditions for this problem are 


1=0, v>0; (6) 
I>0. (7) 


q=P for 


lim g=0 for 

v—0 
The first condition expresses the fact that the 
zeolite at the end of the tube must be saturated, 
since it is in equilibrium with solution containing 
M at a concentration c°. The second condition 
states that there can be no adsorbed material a 
finite distance down the tube if only an infini- 
tesimal volume v has entered the tube. The 
boundary condition for v=0, /=0 is undefined, 
since g changes discontinuously at /=0 when v 
passes through zero. There is an additional 
subsidiary condition 


a) 


f gdl = cv, (8) 
° 


which expresses the conservation of M. 


5 See, for example, E. B. Wilson, Advanced Calculus 
(Ginn and Company, New York, 1912), p. 268. 
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The boundary condition (6) tells us that we 
must take g=P for all negative values of the 
argument in (5). Boundary condition (7) tells us 
that g=0 for all positive values of this argument. 
Placing the argument equal to zero gives an 
explicit equation for the variation of g between 
the limits P and zero. The desired solution is 


therefore 











Kev 
gq=P; I< : (9a) 
P 
1 KPcv\} 
ae ee decd 
1-K l 
Kev cy 
<l<—, (9b) 
KP 
“— (9c) 
q=0; I>—. 9c 
KP 


Figure 1 (calculated for K=0.5) illustrates the 
variation of g with / after a given volume v has 
passed completely through the tube. As K in- 
creases from zero toward one, the curved portion 
becomes more nearly vertical; for K>1, Eq. (9b) 
gives a curve with positive slope, and g becomes a 
triple-valued function of /. Since this would have 
no physical significance, we conclude that g is a 
discontinuous function of / for K2=1. Equations 
(9) are valid for K <1; for K2=1 we have 


g=P; l<cq/P, 
l>cov/P, 


(10a) 
q=0; (10b) 


which gives the dotted curve in Figure 1. The 
solutions (9) and (10) satisfy the conservation 
boundary condition (8). It should be noted that 
the break in the g versus 1 curve comes at the 
point where a fraction K of the total original 
solute has been adsorbed. 

TWO COMPONENTS 


We assume that two substances are being 
adsorbed by equilibrium processes of the type 


HZ + MxX=MiZ + HX, 
(P—qm—@) (4) (qm) (a°+a°-a-a@), 

HZ + M.X¥=MoZ + HX, 
(P—a-—q) (c) (@) (a°+a°—a—a@), (11) 
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where the notation is the same as before. As a 
slight algebraic simplification, we make the 
substitutions 


yi=qi/P; ni=c;/(e1°+c2"); 
x=1P/(c;°+c.°), 


that is, we measure adsorbed material as a frac- 
tion of the total capacity, dissolved material as 
a fraction of the total original dissolved material, 
and length in terms of the length required to 
adsorb completely the dissolved material in unit 
volume of the original solution. Equations (1) 
are then 


On, /dx+0yi1/dv=0; On2/dx+dy2/dv=0, (13) 


and the equilibrium constants for reactions (11) 
are 


(12) 


yi(l —n\—N2) 


my(1—yi—yo). 


2 yo(1—m1—M2) 








2=— . (14) 

n2(1—y1—ye) 
Eliminating the n’s from Eqs. (13) by use of the 
equilibrium equations leads us to the following 
differential equations for the y’s 


Ovi. OY 
[Ke+(1—Ke)y2 -—— (1 —Ke2)y1— 
Ox Ox 

A? dy 
—=0, (15a) 
KiKe ov 





: : de ~. ON 
[Ki+(1 —K,)yi #-- (1 — K1)y2— 
: Ox Ox 


A? Ye 
—=0, (15b) 
Ki Ke Ov 





where 
A=K2K,+K2(1—Ki)yi+- Ki(1—Ko) yo. 


We need solutions of these simultaneous equa- 
tions which will satisfy the boundary conditions 
and, in addition, reduce to the one-component 
equations for the special case K;= Ky». This last 
requirement suggests that we should look for 
special solutions which are similar to those in 
Eqs. (9). As boundary conditions we will use (6) 
and (8), which become, for this problem 








Kyn;° 
m=y'= ; 
Kyn;°+ Kon." 
(16) 
Kon.° 
y2=y92°= 


’ 
Kyn;°+Kon,° 


for x=0, v>0; and 


D 


f ydr=niv; f ysdr=ns%, (17) 
0 0 


c,° C2° 


e _ 
’ Ne aaa . 
€1°+¢2° c1°+c2° 


where 








n= 


The solution y:=constant, ye=constant is obvi- 
ously valid. Inspection of Eqs. (15) shows that a 
change in y; is accompanied by a change in yo, 
and conversely, so that it is not possible to keep 
one y a constant (not including zero) and vary 
the other. If neither y is a constant, a possible 
solution is 


yi° v } 

ae ee 
a 1 e 
2° 

= 


\3 
" (<-) -Ki}, (18b) 
1—Kz, x 


B= Kyn;°+ Kon’. 











where 


There are other special solutions of this form but 
none of them will satisfy the boundary conditions 
and reduce to (9b) for Ki=Ky. The boundary 
conditions show that both y’s remain constant at 
their original values until a value of x is reached 
at which Eqs. (16) and (18) give the same value 
for the y’s; the y’s are given by Eqs. (18) beyond 
this point. The desired solution is thus (for 
Ki <1, K2<1) 


m=", Yo= 2", x < Bu, (19) 


and the values given by Eqs. (18) for «> fv. At 
x = Bv a fraction K,; of component 1 and a fraction 
Kz of component 2 have been adsorbed. Let us 
now assume that K,>K;,. The y’s are given by 
Egs. (18) until the point x = 6v/K,? is reached; at 
this point ye=0 and component 2 has been 
adsorbed completely. We therefore have 


yo=0; x > Bu / 2”, (Ke>K)). (20) 


At this point y2 is continuous but has a discon- 
tinuous derivative. There will therefore be at 
least a discontinuous derivative in the y: curve 
at this point. For x >6v/K.?, y; is determined by 
the differential equation obtained by placing 
y2=0 in Eq. (15a). This gives a one-component 
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Fic, 2. Distribution of adsorbed material: both components 
weakly adsorbed. 








equation 
” Oy, . . Oy1 
K,—+[Kit+(1—K,)y: }}—=0, (21) 
Ox v 
with solutions y;=constant or 
1 v\} 
y= (xe) -Ki}. (22) 
1—K, x 


Boundary condition (17) shows that the y; curve 
becomes horizontal at the point where y2 be- 
comes zero, retains this value until it meets the 
sloping curve (22), and then follows that curve 
down to zero. We have 








Ko—-K, B\? v 
2 ———; pu/K<x<(—) mn, 35) 
1—K, Keo/ Ky, 
aA(*=) -*] 
ni O-) —Kit; 
go OFS ae 
B\2 v v 
( ) —<x<—, (24) 
Ke Ky Ky 
yi=0; x>0/K,. (25) 


At this last point component 1 has been adsorbed 
completely. Figure 2 shows the variation of 4; 
and ye with x after a given volume v has passed 
through the tube. (The actual curves are drawn 
for 2;°=n°=0.5, Ke=0.7, K,=0.5.) 

As Ky increases the curved portion of the y> 
curve becomes more nearly vertical. Just as 
before, if K2=1, the ye curve is discontinuous. 
The boundary condition (17) then shows that the 
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Fic. 3. Distribution of adsorbed material: one component 
strongly adsorbed, one component weakly adsorbed. 


desired solution for the problem K.>1, K,<1 is 











Yo=yo"; x<Bv/Ke, (25) 
y2=0; x>Bv/Ke, (26) 
yi KF > x < Bu ‘Ko, (27) 
Ko—K, B =v 
yi=y1" —-; 0/K.<x<(—) fae (28) 
—K, 2 Ky 
1 7. \ea 
ne tl(nt)'-), 
1—K, x 
B\2 v v 
( ) —<x<—, (29) 
Ke K, Ky 
Vv 
yi=0; x>—. (30) 
Ky 


Equation (28) can be determined from the equi- 
librium conditions at the point x= $v/Keo, al- 
though the boundary condition (17) gives a 
unique solution. Equations (25) to (30) are illus- 
trated by Fig. 3. (The actual curves are drawn 
for n\°=n2°=0.5, Ke=1.5, Ki=0.5.) Equations 
(28) and (29) must be modified if (K2—K,) is 
sufficiently large. The maximum value of 4, is 1, 
therefore, if 








Ko—Ky, 
apa >1, (31) 
1—K, 
Eqs. (28) and (29) must be replaced by 
yi=1; Bv/Ko<x<Ky, (32) 
1 Vv 2 Vv 
1—K, x Ky 
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For ;°=n2°=0.5, condition (31) becomes 
K,>K;,/(2Ki—1), which can be satisfied if 
K,>0.5. This behavior is illustrated in Fig. 4 
(drawn for 2;°=n.°=0.5, K2=4.0, Ki=0.8). 

As K, increases, the last segment of the 4: 
curve approaches the vertical, and for K;=1 we 
have a discontinuity at x =v. For K2>K,21, the 
solutions become 


Yi=1°, Yo=2"; eee, 
Ke 
Bu (34) 
yi=1, y2=0; a SEP, 
Ko 


yi=0, y2=0; x>v. 


MORE THAN TWO COMPONENTS 


If we have m equilibria of the type (11), with 
equilibrium constants defined by 


m 


yi(1— 2 ni) 
i=1 
kK ,;=———_ ; ¢=1,2,---m, (35) 
n(1—> yi) 
i=1 





the set of partial differential equations becomes 


dy; oA Kid? dy; 
Se ce 
Ox Ox Kiko: . -Km Ov 
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where 


of a flake 
a= (KiKs--Kx)| 1+ (- )a} 


i=1 





with the boundary conditions 


Kn; 4) 
ye =- a f y dx =n ;"v, 
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m 


ye K in e 
i=1 





7=1,2,---m. (37) 
Solutions of Eqs. (36) of the form (18) are 


yi° m v b 
Vi= iT % Kn?)-] -Ki. (38) 
1-K; i=1 x 


If all the K’s are less than 1, the y’s will remain 
constant at the values y;° until the point 





™m 


x=(>> Kin;)v 
i=1 


is reached. At this point a fraction K; of each 7’th 
component has been adsorbed. The y’s would 
then follow Eqs. (38), giving curves similar to 
the first part of Fig. 2 but with m curves instead 
of two, until the component with the largest K 
has been adsorbed completely. The problem is 
then reduced to one of lower order and the same 
general type of variation is repeated, giving a 
succession of curves similar to those of Fig. 2. If 
some of the K’s are greater than unity, the 
corresponding y’s will be discontinuous, and we 
will have saturation of the adsorbing material up 
to the point at which all components with K’s 
greater than unity have been completely 


adsorbed. 


CONCLUSION 


The above results apply to the original deposi- 
tion of the adsorbed material. In the process of 
chromatography, the chromatogram is ‘“de- 
veloped” by pouring pure solvent down the tube. 
The region containing the adsorbed material is 
then displaced along the tube. For our particular 
problem, this process corresponds to pouring HX 
at concentration (c;°+c2°) down the tube. The 
same equations apply to this problem, but they 
apply to both the front and rear boundaries of the 
adsorbed material. Equation (9b) has a positive 
slope for K>1, and therefore is meaningless 
when applied to the front boundary, but not 
when applied to the rear boundary. Therefore, 
for K>1, the front boundary is sharp during 
development but the rear boundary is not; for 
K <1 the rear boundary is sharp and the front 
boundary is not.”* In developing a pattern such 
as that of Fig. 3 (K.>1, Ki<1), the y: curve 
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Fic. 4. Distribution of adsorbed material: both components 
strongly adsorbed. 
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Fic. 5. Concentration of adsorbed material as a function 
of concentration of solute, illustrating type of curvature of 
the isotherms. 


will have a sharp rear boundary, and the y2 curve 
will develop a diffuse rear boundary. The region 
over which the two curves overlap shortens and 
eventually disappears. For the particular situa- 
tion illustrated in Fig. 3, the solution during 
development which reaches the point at which y»2 
becomes zero will have the same concentration as 
had the original solution, and this point will 
therefore move a distance BV /K2, where V is the 
volume of solvent used in development. The rear 
boundary of yi, which is sharp, will move a 
distance BV /K,, since a volume V of solvent will 
remove 7,°V units of component 1 from a region 
where the concentration of adsorbed material is 
yi°= K,n,°/8. The two components will therefore 
be separated completely when 


a 


Ki 
V= ee ee (39) 


27441 


B(V+)/Ke=BV/K,, or 


where vp is the volume of solution used in the 
original adsorption. If both K’s are greater than 
one or both less than one, complete separation 
will be more difficult. 

The explicit formulas for the amount of ma- 
terial adsorbed as a function of distance down the 
tube are of some interest. The total adsorbed 
material in a length L of the tube is 


L 
m= f qidl. 
0 


(40) 
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Let L=alpo, where Lo is the length required to 
adsorb all the solute in a volume vp of solution. 
For the one component problem, with K <1, we 
have Lyo=c%/KP. The fraction of the original 
solute adsorbed in a length L is then 


Q 





—s ask’, (41) 
K 
2a!'—a—K 
=< -; K*Sasl1. (42) 
1-K 


Although the above algebraic relations apply 
only to a special type of adsorption, the results 
should have general qualitative validity. If 
Eq. (14) is solved for y:;, we have 


yi=Kyn,/[1—(1—Ki)m1—(1—Ke) ny]. (43) 


The behavior of this function for the region of 
small ; (holding m2 constant) is shown in Fig. 5. 
K,=1 gives the dividing line between substances 
which give a sharp front boundary and those 
which give a diffuse front boundary. As shown 
above, this criterion is independent of Ks. The 
important feature is not the actual slope of the 
y1 versus nN; Curve, since this depends on (1 — K2)n2 
as well as on K,, but the way in which the slope 
changes with increasing m;. If the adsorption 
isotherm is represented by a more general equa- 
tion than Eq. (43) we expect a diffuse front 
boundary .if the slope of the y,; versus n, curve 
increases with increasing m, for small 2, and a 
sharp front boundary if the slope of the curve 
decreases with increasing ”; in this region. This 
conclusion is in accord with the original sugges- 
tion of DeVault.? Of course, the exact shape of 
the diffuse boundary will depend on the particular 
adsorption isotherm. 

This problem was suggested to the author by 
Dr. F. R. Duke. The author wishes to acknowl- 
edge his indebtedness to Dr. Duke for the 
formulation of the problem and for valuable 
discussions regarding the validity of the as- 
sumptions made and the results obtained. 
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HE present paper is a survey of the appear- 
ance in the electron microscope of a number 
of important catalysts. Its purpose is to delineate 
the forms of catalytic materials and to serve as an 
introduction .to further quantitative studies. 
Previous work on the subject has been reported 
in Germany by M. V. Ardenne and D. Beischer.' 
The electron micrographs were taken on the 
RCA B Electron Microscope at the RCA Re- 
search Laboratories at Princeton. The catalyst 
material was mounted for examination in the 
following way: The catalyst was triturated with 
a two percent solution of pyroxylin in amyl 
acetate. A drop of the suspension so obtained was 
placed on the clean surface of distilled water and 
allowed to spread freely to give a thin film about 
200A thick. Two hundred mesh supporting nickel 
screens were then placed on the collodion, the 
film and screen picked up, dried, and examined in 
the microscope. The grinding with the pyroxylin 
solution seems to break up the secondary aggre- 
gates or clumps, and the spreading of the film on 
water often facilitates the separation of one 
particle from another. 
The results obtained are presented in a series 
of electron micrographs. The scale is indicated on 





Fic. 1. Platinum oxide. 


1M. V. Ardenne and D. Beischer, Angew. Chem. 53, 103 


(1940). 


the side and it must be recalled that 1 micron is 
equal to 10,000A. Figure 1 represents platinum 
oxide prepared according to the directions of 
Adams.” It can be seen that there are flakes of 
material ranging in size from a tenth of a micron 
to two microns. In addition to these more obvious 
particles one should note that the whole back- 
ground of the collodion membrane is filled with 
finely-divided, homogeneous material whose par- 
ticle size is near the limit of the resolution of this 
micrograph (about 200A). This finely-divided 
material is interpreted to be the platinum oxide 
from which the platinum catalyst is to be made 
by reduction. This interpretation is to be con- 
firmed by subsequent electron diffraction study. 

Figure 2 represents the platinum catalyst ob- 
tained by the reduction of the previously de- 
scribed platinum oxide. The dark masses are the 
platinum particles. It is seen that they range in 
size from 200 to 10,000A. There is definite 
evidence from the straight sides on some of the 
platinum particles that they are crystalline. The 
micrograph shows several folds and holes in the 
supporting collodion membrane. 

The next series of micrographs deals with the 
appearance of asbestos and platinized asbestos. 





Fic. 2. Adams’ platinum catalyst. 


2 Gilman, Organic Syntheses (John Wiley & Sons. Inc., 
New York, 1932), p. 53. 
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Fic. 3. Asbestos fibers. 


Figure 3 represents a dispersion of asbestos 
fibers. The breadth of these fibers is 3000 and 
6000A and their straight edge is particularly 
noteworthy. The ovals in the picture are merely 
holes in the supporting collodion membrane. 
Impregnation with chlorplatinic acid does not 
change the appearance of the fiber—the chlorpla- 
tinic acid seems to be very finely dispersed in the 
fiber. 

Figure 4 presents the appearance of an active 
platinized asbestos catalyst. This particular prep- 
aration has had considerable use in catalyzing the 
hydrogen-oxygen reaction. The formerly clear 
asbestos fiber is thickly covered with material 
which is undoubtedly finely-divided platinum. A 
significant feature of this micrograph is the 
cubical crystal (one micron on side) of platinum 
which is attached to the fiber by one side and 
corroded at the two adjacent sides. Other micro- 
graphs indicate that the cubical shape and this 
type of attachment to the fiber, is typical of the 
platinum. Figure 5 is another electron micro- 
graph of a platinized asbestos catalyst. The field 
here chosen represents very thin slivers of asbestos 
which are partially transparent to the electrons. 
The dense aggregates disclose the size and distri- 
bution of the platinum material on this support. 
One can detect in this micrograph particles very 
close to 50A in diameter, the very resolution 
limit of the microscope. 

Figure 6 shows the appearance of the catalyst 
after it has been subjected to a sintering process 
of heating to 600°C. The appearance of this 


Fic. 5. Platinized asbestos. 


inactive catalyst is in striking contrast to the 
previous platinum photographs. The edge of the 
asbestos is clean, the fine particles of platinum 
have disappeared and the main bulk of the 
catalyst has formed a swarm-like aggregate on 
the fiber. There is also some indication that the 
corners and edges of the material have rounded 
off. 

Figure 7 represents ten-percent platinum on 
activated charcoal. The characteristic feature of 
this catalyst is the plate-like structure of the 
support which is particularly evident under and 
beside the large piece of charcoal. The plate form 
of charcoal is consistent with its large surface 
area of about 1000A per gram.’ The platinum is 


3P, Emmett, in Advances of Colloid Science (Interscience 
Publishing Corporation, New York, .1940) 1, 29. 
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so well dispersed that in spite of the fact that 
there is two percent of it on the charcoal, one can 
see the dense platinum aggregates only here and 
there on the field. 

Figure 8 shows the appearance of Alorco 
activated alumina. It is seen that the product is 
homogeneous and consists of some particles of 
100A and less. The surface of this material as 
determined by the Emmett-Brunauer gas ad- 
sorption method is equal to 270 sq. meters per 
gram. Steaming at 600°C for 40 hours reduces the 
surface to 13 sq. meters per gram. This product 


Fic. 6. Sintered platinized asbestos. 


Fic. 7. 10-percent platinum on charcoal. 


. 
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Fic. 9. Steamed alumina. 


Fic. 10. Silica gel. 
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Fic. 11. Alumina silicate cracking catalyst. 


Fic. 12. Zinc oxide methanol catalyst. 


is shown in Fig. 9. The particles are considerably 
larger and have a definite plate like structure. 

Active silica gel is shown in Fig. 10. In addition 
to the larger clumps the field is full of very small 
particles whose diameter is of the order of 200A. 
This is consistent with the large specific area of 
the silica gel. 

The electron micrograph of an alumina silicate 
catalyst used for catalytic cracking of petroleum 
is shown in Fig. 11. There is marked evidence for 
spheres of 200A-diameter and rod-like particles 
1000A long and 300A wide. 

The appearance of zinc oxide catalyst made by 
the decomposition of zinc oxalate and used for 
methanol synthesis is shown in Fig. 12. The 
particles are very finely divided having diameters 
which fall in the range of 100-200A. Their 


Fic. 13. Fischer synthin catalyst. 


spherical shape may be due to the inability of 
the microscope to resolve the shape in this par- 
ticular micrograph. The appearance of this zinc 
oxide catalyst is in marked contrast to the zinc 
oxide smoke which has the well-known needle like 
structure. 

The last electron micrograph (Fig. 13) to be 
presented is that of the Fischer catalyst con- 
sisting of nickel-thoria supported on kieselguhr. 
This catalyst is used for the formation of pe- 
troleum hydrocarbons from carbon monoxide and 
hydrogen. The kieselguhr is seen to be composed 
of diatoms which are particles full of holes 1300A 
in diameter. These openings offer a large surface 
for the support of the nickel and thoria. It can be 
seen from the picture that some of these openings 
are completely filled while others are only 
partially filled. 

The electron micrographs presented give a 
survey of the various ways large surfaces are 
produced by the art of the catalytic chemist. One 
method consists of using finely divided materials 
as in the case of active alumina, silica gel, etc.; 
the second available method is to take advantage 
of the rod-like nature of asbestos fiber as a 
support for metallic particles like platinum. 
Another method consists in using as a support 
the plate-like structure of materials such as 
active charcoal. Another way of attaining a large 
surface per unit weight is the use as a support of 
such porous material as the diatoms found in 
kieselguhr. Since catalytic activity is associated 
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with large surface area, the electron microscope is 
a powerful tool for the evaluation of a given 
material as a possible catalyst. 

This investigation has been carried out in part 
at the RCA Research Laboratories at Princeton. 
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A few selected topics from four papers by Samuel on the theory of valence are critically 
examined. The conclusions are drawn that Samuel's bond-moments for molecules containing 
*‘pentacovalent nitrogen” are not as satisfactory as claimed, and provide no evidence in favor 
of the existence of pentacovalent nitrogen; that Samuel’s eigenfunctions for structures with 
pentacovalent nitrogen are identical with the ones set up in the now conventional manner for 
structures without pentacovalent nitrogen but with semipolar double bonds instead; that the 
existence of antibonding orbitals seems to be theoretically necessary; and that several of 
Samuel’s further conclusions are either based upon the use of words in unfamiliar senses, or else 


are erroneous. 


N a recent series of papers by the same name 
as the present one, Samuel! has come to con- 
clusions which appear to be at variance with the 
present generally accepted theories of valence. 
He proposes that many of the features of modern 
theory be abandoned, and that a return be made 
to the corresponding features of classical theory.* 
For example, he considers that the semipolar 
double bond and the octet rule are without either 
experimental or theoretical proof and, indeed, 
that their acceptance makes more difficult the 
interpretation of a number of the physical prop- 
erties of the substances concerned. In order to 
support his thesis, he advances two types of 
argument. First, on the basis of the classical 
structures, he establishes sets of constants, with 
the use of which, in his opinion, the molecular 


1R. Samuel, J. Chem. Phys. 12, 167, 180, 380, 521 (1944). 

*For lack of better terms, the words “‘classical’’ and 
“modern”’ will be used throughout this paper to describe, 
respectively, the theories, structures, etc. which Samuel 
adopts, and the corresponding now generally accepted 
theories, structures, etc. which he rejects. The distinction 
between the classical and the modern points of view will 
become more apparent as we proceed. In the use of this 
terminology, no judgment regarding the relative merits of 
the two viewpoints is intended. 


parachors, dipole moments, and the like can be 
calculated more satisfactorily than with the use 
of any other sets of analogous constants that have 
been derived hitherto on the basis of modern 
theory. And second, he points out several respects 
in which, in his opinion, the quaantum-mechanical 
arguments that have been used to ‘‘prove”’ the 
modern theory are inconclusive, illogical, and 
contrary to spectroscopic evidence. 

Since the significance of Samuel’s papers may 
be somewhat obscure to the reader who is not 
particularly well acquainted with the field, | 
wish here to discuss critically, and in some detail, 
a selected few of the salient points which Samuel 
makes. 


A 


As has already been mentioned, Samuel is con- 
cerned to a considerable extent with the use of 
parachors, dipole moments, and the like, in the 
study of the finer details of molecular structure. 
He points out, quite correctly, that it is impossi- 
ble, on the basis of such properties, to prove the 
existence of semipolar bonds, which he does not 
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accept, to justify the octet rule, or to support any 
of several other features of modern valence 
theory. I wish here, however, to emphasize more 
strongly than Samuel does, that it is equally 
impossible, on the basis of these same properties, 
to establish the pentacovalence of nitrogen, which 
Samuel adopts, to disprove the octet rule, or to 
support the classical structures. The most that 
one can hope to do with the use of empirically 
obtained atomic parachors, bond-moments, and 
the like, is to make it seem probable that a given 
bond in a given molecule is (or is not, as the case 
may be) of the same type as a second bond in a 
second molecule. Whether the bond in question 
is double, semipolar double, hybrid, or whatnot 
cannot be decided in such a way, but only on the 
basis of theoretical arguments and with the use 
of a self-consistent set of definitions. This point 
will be elaborated further in the present section 
and in the following one. 

Let us begin the discussion by considering 
Samuel’s calculation of the dipole moments of 
several types of compound containing ‘‘penta- 
covalent nitrogen.’’ With the assumption of the 
correctness of the classical structures, in which 
the covalence of nitrogen is permitted to be either 
three or five, as required, Samuel derives a set of 
bond-moments, which he considers to be constant 
from molecule to molecule. In tables (and also in 
the text and in diagrams illustrating the vector 
summation of bond-moments), he lists in his 
third paper a number of such values, of which the 
following ones will be important for our present 
purposes. 


H—C,, 0.3D 
H—C,, 0.4D 
Cy—N™ 0.5 D 
C..—N"™' 0.9 D 


NY=O0 44D NY=C 2.6D 
NY—O 14D NY—C 1.0D 
NY=SN'™"4.2D Cyur—Cy, 0.6 D 
NY=N"! 2.8 D 


In each case, the superscript designates the 
covalence of the nitrogen atom, and the sub- 
scripts “al” and ‘‘ar’’ mean ‘aliphatic’ and 
“aromatic,” respectively. The atom given at the 
left of each bond is considered to be at the 
positive end of the corresponding dipole. These 
values are apparently intended to be used in 
conjunction with the dipole moments measured 
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in solution, rather than in the gas phase, since the 
gas-phase moments are not known for most of the 
compounds considered. 

As Samuel shows, the above bond-moments 
permit the calculation of the molecular moments 
of a considerable number of different types of 
substance. The agreement between calculated 
and observed values, however, is not always as 
satisfactory as the reader might infer from 
Samuel’s discussion. The following examples will 
illustrate the magnitudes of the discrepancies in 
the worst cases. The most serious of these dis- 
crepancies occurs with trimethylamine oxide, for 
which the calculated and observed? moments are 
3.7 D and 5.0 D, respectively. Samuel does not 
discuss the moment of this substance. In some 
cases, the satisfactory agreement which Samuel 
obtains between calculated and observed mo- 
ments appears to be due to numerical errors. 
For example, Samuel calculates the dipole mo- 
ment of azoxybenzene to be 1.9 D, which is only 
0.2 D greater than the observed value of 1.7 D. 
However, if the bond-moments and bond-angles 
are those given explicitly in his vector diagram, 
the resultant calculated moment is actually 2.3 D. 
The error is therefore three times as great as 
stated by Samuel. Furthermore, the calculated 
and observed moments of phenyl azide are 2.3 D 
and 1.55 D, respectively. For this substance, 
Samuel makes no calculation ; instead, he merely 
states, without reference to specific data, that the 
moment of the azide group proper is 1.4 D to 
1.6 D, as compared with the calculated value of 
1.4 D. (Actually, in footnote 19 of his second 
paper,! Samuel refers to a paper*® in which 
Sidgwick states that the moment of the azide 
group proper is ‘‘ca. 0.”) And finally, Samuel 
gives 6.4 D and 1.4 D as the calculated moments 
of the a- and 8-N-methyl-p-nitrobenzophenone 
oximes, respectively, as compared with the ob- 
served values of 6.6 D and 1.2 D, respectively. 
However, the correct calculated values,* obtained 


2E. P. Linton, J. Am. Chem. Soc. 62, 1945 (1940). 

3N. V. Sidgwick, Chem. Rev. 19, 183 (1936). 

* According to a statement received from Dr. Samuel, 
the value of 132° given for the N—C—C bond angle in his 
diagram on page 387 is a misprint for 122°. Similarly, the 
value of 0.8 D for the moment of the bond between the nitro 
N and the phenyl C in the same diagram appears to be a 
misprint for 1.0 D. These corrected values are used in the 
calculations referred to above. 
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with the use of his bond-moments and bond- 
angles, are instead 6.0 D and 1.7 D, respectively. 

Let us now return to Samuel’s list of bond- 
moments. If we consider that all the moments 
which do not involve the ‘‘pentacovalent nitrogen 
atom,’’ NY, are correct, we are left with altogether 
six independent parameters, namely, the mo- 
ments of the NYO, NY—O, NY=N"!, NY==N!1, 
NY=C, and NY—C bonds. Consequently, Samuel 
should be able to calculate with no error the 
moments of six independent types of molecule, 
regardless of whether the bond-moments are even 
approximately constant or not. Actually, Samuel 
considers altogether nine types of molecule, 
namely, nitro compounds, methy] nitrate, nitrous 
oxide, azides, aliphatic diazo compounds, azoxy- 
benzenes, iso-azoxybenzenes, and a- and B-N- 
methyl oximes; and he ignores a tenth type, 
namely trimethylamine oxide. It may be ques- 
tioned whether the azoxybenzenes and the iso- 
azoxybenzenes are really independent of each 
other, since, being stereoisomers, they contain 
exactly the same kinds of bond and differ only in 
bond-angles. Similarly, and for the same reason, 
it may be questioned also whether the a- and the 
8-N-methyl oximes are independent. In any 
event, however, Samuel’s method of calculation 
can be applied, at the very most, to only ten 
different types of molecule containing ‘‘penta- 
covalent nitrogen.’’ Of the ten corresponding 
calculations, those for phenyl azide and for 
trimethylamine oxide are in error by 0.75 D and 
1.3 D, respectively, as we have seen, and those 
for azoxybenzene and for a- and 8-N-methyl-p- 
nitrobenzophenone oxime are, at best, not very 
good. It is evident, therefore, that Samuel’s 
calculations are not sufficient to establish that 
the bond-moments which he has derived are 
more than very roughly constant from molecule 
to molecule. It is even more evident that these 
calculations provide no support for Samuel’s 
belief that the bonds in question involve penta- 
covalent nitrogen. 

Before leaving the present discussion, let us 
consider the modern interpretation of the above 
substances. Samuel’s NY==O bond, for example, 
occurs in the nitro compounds, nitrous oxide, the 
azoxybenzenes, the N-methyl oximes, methyl 
nitrate, and trimethylamine oxide. For these 
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molecules, Samuel writes the classical structures 
O—N=—O, C,H;—N=N " 
| N=N=—O, | 
R CH; 


R CH;—O CHs 
| 


| | 
C=N—0, N=0, and CH;—N=0, 


ta | | 
R CH; O CHs 


respectively. On the other hand, from the modern 
point of view, the structures must be written 
instead as 
-O—N*=0 O=Nt—O- 
| “* | 


| 
R R 
-N=Nt=—0eN=N*t—O- 


N- “a =f) C.H; N=Nt—O- 
“> 


| | 
CeHs C.Hs 


R R 
-¢_N+040C=N—0- 
k CH oR CH, 

CH;—O 


CH;—O+ CH;—O 


| | | 
N+*=0— N+—O-@— N+—O- 
| | | 


ez 


| 


O- O 


CH; 


| 
CH;—Nt—O- 
| 


| 
CH; 


respectively.* The double-headed arrows here, as 
usual indicate that resonance occurs between the 
structures stated. In all of the first four molecules, 


*Here, and throughout the remainder of the present 
paper, a semipolar double bond between two atoms A and B 
is represented as At—B~-, and not as AB. The former 
mode of representation, which appears now to be quite 
generally adopted in this country, has the double advantage 
of being unambiguous and of suggesting correctly that the 
semipolar bond At—B- is an ordinary single covalent bond 
A—B plus an ionic bond At+B~; on the other hand, the 
latter mode of representation, which is employed by 
Samuel, could be confusing since it suggests erroneously 
that the bond A=B is a single bond A—B plus a singlet 
linkage A—B. 
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in which the resonance is between just two struc- 
tures, the nitrogen-oxygen bond under discussion 
is seen to be of more or less similar, but of course 
not identical, type. Consequently, we can hardly 
be surprised if the moment of this bond proves 
to be sufficiently nearly constant to make possible 
the rather limited success of Samuel’s method of 
calculation, since, as has already been seen, fairly 
large deviations from constancy would apparently 
not be inconsistent with the experimental data. 
In methyl nitrate and trimethylamine oxide, 
however, the situation is different, and the 
“‘NY—O bonds” are not very closely analogous to 
those in the other molecules, since, in the former 
substance, the resonance includes three structures* 
and, in the latter, no resonance is presumed to 
occur at all. Consequently, the modern viewpoint 
might lead us to expect that here, if anywhere, 
relatively serious discrepancies between calcu- 
lated and observed dipole moments would be 
found. With methyl nitrate, this expectation can 
be neither verified nor shown to be incorrect, 
since this substance is the only one considered by 
Samuel to contain the NY—O single bond ; conse- 
quently, perfect agreement between the calcu- 
lated and observed moments can be obtained by 
suitable adjustment of this NYN—O bond-moment, 
no matter what values are assigned to the re- 
maining bond-moments. With trimethylamine 
oxide, on the other hand, the expected dis- 
crepancy does indeed occur, as has already been 
noted. 

In the present connection, mention may be 
made of the fact that the relative contributions 
of resonating structures are riot as sensitive to the 
differences between their energies as is often sup- 
posed by people who are unfamiliar with the 
subject. In particular, no analogy exists with the 
way in which the position of equilibrium in a 
tautomeric system depends upon the difference in 
free energy between the isomers. For example, 
calculation shows that, if a molecule resonates 
between two structures that differ in energy by 
even as much as 100 kcal. per mole or more, there 
is no reason why the less stable structure might 
not make, say, half as great a contribution as 

* Actually, methyl nitrate may not be so different 
from the preceding substances after all, since the second 
of its above structures is doubtless relatively unstable on 


account of the unfavorable charge distribution. This fact 
does not affect the following discussion, however. 


does the more stable one. This is because the 
relative contributions of the structures are de- 
termined not only by the difference between their 
energies but also by the magnitude of their 
interaction integral (i.e., of the non-diagonal 
matrix element of the energy). In striking con- 
trast, however, a similar difference of 100 kcal. 
per mole between the free energies of two 
tautomers would necessarily correspond to an 
equilibrium constant of about 10”, so that at 
equilibrium the less stable form could be said not 
to be present at all. In his fourth paper,! in which 
he criticizes the modern interpretation of the 
moments of the above-mentioned substances, 
Samuel seems to overestimate both the im- 
portance of the variations in the relative contri- 
butions of the resonating structures, and also the 
success of his method of calculation. 


B 


In this section we shall, like Samuel, consider 
whether the fundamental quantum-mechanical 
theory of valence permits a decision to be made 
between the classical and the modern structures. 
We shall find it necessary first, however, to 
examine rather carefully the exact meanings of 
the terms and symbols used. 

Let us begin by considering what quantum- 
mechanical meaning is to be attributed to the 
statement that trimethylamine oxide, for ex- 
ample, has the classical structure I and not the 
modern one II (or vice versa). Obviously, it is not 
imagined that with a sufficiently powerful “‘y-ray 





CH; CH; 
| | 
CH;—N=O CH;—Nt+—O- 
| | 
CH; CH; 
I II 


microscope”’ one would actually be able to see the 
bonds, either as straight lines or as pairs of dots 
in the electronic representation. Neither, obvi- 
ously, is it imagined that the distinction could 
be made on the basis of the chemical properties 
of the substance; actually, these would be even 
less useful for this purpose than the physical 
properties considered in section A. The view, 
which is probably shared by the majority of the 
people who have given recent critical thought to 
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the question, seems to be approximately the one 
developed in the remainder of this section. 

The eigenfunction of a molecule, which defines 
the distribution of charge within it, provides the 
most precise and the most nearly complete 
description of the molecule permitted by quantum 
mechanics. Unfortunately, however, this eigen- 
function cannot be determined exactly for such a 
complicated molecule, and, even if some well- 
disposed deity were kind enough to tell us what 
it is, the original problem would still be unsolved. 
The reason for this is that the correct eigen- 
function would contain no recognizable feature 
which could be identified either as a double bond 
or as a semipolar double bond between the 
nitrogen and oxygen atoms. It would give the 
distribution of charge and, as corollaries, the 
geometry of the molecule, the magnitudes and 
directions of the restoring forces opposing the 
displacements of the atoms from their most 
stable positions, and so on, but it would say 
nothing about the types of bond present. 

In order to obtain information about the types 
of bond, one must leave the rigorous quantum- 
mechanical procedure and must introduce approx- 
imations. Two different approximate methods 
have been adopted for effecting a correlation be- 
tween quantum mechanics, on the one hand, and 
structure in the usual chemical sense, on the 
other. The first of these is known as the valence- 
bond, or pair-bond, or Heitler-London-Slater- 
Pauling (HLSP) method, and the second as 
the molecular-orbital or Hund-Mulliken (HM) 
method. In either method, a procedure exists by 
which an approximate eigenfunction can be set 
up to correspond to any given structure. If there 
is adequate reason to believe that an eigenfunc- 
tion set up in either of these ways is a fair ap- 
proximation to the correct, but unknown, eigen- 
function of a given molecule, then, by definition, 
the molecule is said to have the corresponding 
structure. On the other hand, if there is reason to 
believe that a linear combination of two or more 
such functions should give a materially better 
approximation to the correct eigenfunction than 
any one of them alone, then, again by definition, 
the molecule is said to be a resonance hybrid of 
the two or more corresponding structures. It 
should be clear that, since no single approximate 
function which is set up in accordance with either 





of the two methods can be exactly correct, every 
molecule is, strictly, a resonance hybrid of a 
number of structures. In fact, it can be shown 
that a perfect representation of even the simplest 
molecule could be obtained only by consideration 
of resonance among all the structures that can be 
imagined, including an infinite number of them 
that cannot be represented by any of the con- 
ventional chemical symbols. However, it is a 
fundamental assumption of both the pair-bond 
and molecular-orbital treatments that any mole- 
cule can be described, with a precision that is 
sufficient for most chemical and spectroscopic 
purposes, by means of a single structure in the 
above sense, or, at worst, as a resonance hybrid 
of a small number of such structures. (In the 
interests of strict accuracy, the following quali- 
fications should be noted. Some structures corre- 
spond to eigenfunctions which are identically 
equal to zero, as will be mentioned again later; 
and some turn out to be merely resonance hybrids 
of two or more other ones. Obviously, structures 
of the first type may always be ignored; and 
those of the second type may also be ignored if 
their various components have already been 
included in the treatment.) 

From what has been said, it follows that the 
answer to the question whether trimethylamine 
oxide has the structure I or II, above, depends to 
a considerable extent upon the way in which 
structure and approximate eigenfunction are 
correlated. It is quite possible that two different 
people, approaching the problem from different 
angles, might relate a given function to two 
entirely different structures. As we shall see, this 
fact seems to be a major reason for the disagree- 
ment between Samuel’s views and the modern 
ones, since Samuel essentially proposes a new 
way of correlating approximate eigenfunction 
and structure. 

Let us now consider for a moment the way in 
which the pair-bond and molecular-orbital treat- 
ments deal with the structure of trimethylamine 
oxide. By either method, one can set up an ap- 
proximate eigenfunction for each of the structures 
I and II, above. The functions for the two 
structures differ in the important respect that, 
whereas the ones for structure II make use of only 
K- and L-shell atomic orbitals, the ones for 
structure I require that also the M-shell orbitals 








244 G. W. 


of the nitrogen atom be taken into account. Since 
the orbitals of the M-shell are much less stable 
than those of the K- and L-shells, the modern 
view is that structure I should be considerably 
less stable than II. Consequently, it is supposed 
that, although the resonance must to some extent, 
of course, include both structures, in addition to 
many others, nevertheless structure II gives 
much the best description of the molecule that 
any single structure can give. 

In this connection, it will be instructive to con- 
sider the reason why the structure I, from either 
the pair-bond or molecular-orbital point of view, 
requires use of the M-shell orbitals on the nitro- 
gen atom. By following the formal rules for 
setting up approximate eigenfunctions to repre- 
sent specified structures, one can indeed write 
down an expression to represent I, in which only 
K- and L-shell orbitals are used. However, this 
expression, after it has been made antisymmetric 
as required by the Pauli exclusion principle, i 
expressible as a sum of terms, of which a number 
are identically equal to zero. If these vanishing 
terms are omitted from further consideration, the 
function which remains proves to be one corre- 
sponding to some modern structure, or to reso- 
nance among several modern structures, without 
any pentacovalent nitrogen atom. (Cf. the dis- 
cussion, below, of the configurations of the 
excited nitric oxide and nitrogen molecules.) 
It is evident, therefore, that, with the restric- 
tion that only K- and L-shell orbitals be used, 
and within the limitations imposed by the defini- 
tions inherent in the pair-bond and molecular- 
orbital methods, the qctet rule is not an added 
arbitrary assumption, but is a necessary conse- 
quence of quantum mechanics and of the exclu- 
sion principle. It cannot be violated even if the 
attempt is made to do so.* 

Samuel, however, considers that structure I, 
with pentacovalent nitrogen, can be obtained 
even without use of the M-shell orbitals, and he 


* With phosphorus, sulfur, and any of the other atoms of 
the second row of the periodic table, the pair-bond and 
molecular-orbital treatments permit violations of the octet 
rule without the use of the orbitals of an outer valence shell. 
This is because, for these atoms, the valence shell, being the 
M-shell, contains 3d, as well as 3s and 3, orbitals. Even 
though the 3d orbitals are relatively unstable, it is, in fact, 
decidedly questionable whether the octet rule is rigorously 
valid, even from the modern point of view, with such 
atoms. The situation is similar with the still heavier atoms 
of the third and subsequent rows. 
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shows how a non-vanishing approximate eigen- 
function can be set up to represent that structure. 
We shall consider, first, the relation between his 
function and the ones which would be set up on 
the basis of the pair-bond and molecular orbital 
methods to represent the modern structure II; 
and, second, the significance of the statement 
that his function actually does represent the 
structure I. 

Samuel explains, but does not write down 
explicitly, the configuration of an excited mole- 
cule of nitric oxide, which contains three unpaired 
electrons and so may be considered to be trivalent 
and capable of forming three covalent bonds to, 
for example, the three methyl carbon atoms in 
trimethylamine oxide. From his discussion,‘ and 
also from analogy with his configuration of an 
excited, divalent nitrogen molecule,' it is evident 
that the configuration under discussion can be 
expressed in his symbols as 


o°(2s), o*(2s), w(2p), o°(2p), x*?*(2p), 


where reference to the K-shell electrons is omitted 
as not being necessary for our purposes. (Actu- 
ally, it makes no difference exactly what con- 
figuration is assumed here, since essentially the 
same discussion as the following could be based 
upom any reasonable configuration.) The three 
unpaired electrons include the one occupying the 
o*(2s) orbital and the two occupying the two 
different 2*(2p) orbitals. Since these electrons 
are to be used to form bonds to the nitrogen atom, 
the orbitals which they occupy must be, re- 
spectively, essentially 2s and 2pz atomic orbitals 
localized about the nitrogen atom ; otherwise, the 
resulting bonds would involve the oxygen as well 
as the nitrogen or, in other words, would be 
resonance hybrids of nitrogen-carbon and oxygen- 
carbon bonds. Consequently, since Samuel re- 
stricts the consideration to the atomic orbitals of 
the K- and L-shells, we may replace the symbols 
o*(2s) and 2*(2p) by 2sy and 2pzn, respectively. 
Moreover, it can be shown by application of the 
Pauli exclusion principle that the localization of 
the orbitals o*(2s) and 2*(2p) about the nitrogen 
atom requires the simultaneous localization of the 
further orbitals ¢(2s) and 2(2p) about the oxygen 


4 Page 385 of reference 1. 
5 Page 184 of reference 1. 
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atom.* These last orbitals can therefore be repre- 
sented as 2so and 2pz0, respectively. Conse- 
quently, the above configuration assumes the 
form 

250", 2sn, 2p70', o*(2p), 2prn’, 


in which only the orbital o(2)) is still a molecular 
one. 

Now, Samuel considers that his configuration 
of the excited nitric oxide molecule corresponds 
to the structure II], so that, when bonds to three 
methyl groups are formed, the structure of the 
resulting amine oxide is indeed the classical one, 
]. However, inspection of the configuration in the 


| | 
—N=0 —Nt-o- 

| | 

Il IV 


second of its two above forms shows that it is 
exactly the configuration which would be set up, 
by a familiar method embodying some features 
of both the pair-bond and molecular-orbital 
treatments, to represent structure IV. On this 
basis, then, the resulting structure of the amine 
oxide is the modern one IJ. Consequently, it 


follows that, when Samuel shows how an eigen- 
function can be set up to represent the classical 
structure I, he is not really contradicting the 
modern ideas of valence, nor is he adding any- 
thing essentially new to them; he is instead re- 
expressing the same ideas in different language 
and with a different symbol. 

With nitrous oxide, and in general with the 
other molecules which, from the modern point of 
view, receive important contributions from two 
resonating structures, the situation is somewhat 
more complicated. In his third paper,’ Samuel 
describes nitrous oxide as being formed from 
NO+N, that is, from the excited trivalent nitric 
oxide molecule discussed above, plus a normal 
nitrogen atom with its three unpaired electrons. 
If the nitric oxide is assigned Samuel’s structure 
III, then the resulting nitrous oxide must be 
represented by the classical structure V. On the 

* The proof of this statement, which is too long to be 
given here in detail, is based upon the fact that all the 
different molecular orbitals in the molecule must be 
mutually orthogonal; if they are not all orthogonal, the 
resulting antisymmetric eigenfunction for the molecule as a 
whole can be broken up into a sum of terms, of which one 


does involve only mutually orthogonal orbitals, and all the 
others are identically equal to zero. 
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other hand, if the nitric oxide is assigned the 
structure IV instead, then the resulting nitrous 
oxide must be represented by the modern 
structure VJ. This last structure is, of course, 


N=N=O N=N+—O- -N=N+=0 
V VI Vil 


only one of the two structures between which the 
modern theory presumes resonance to occur. The 
second modern structure VII arises automati- 
cally, however, in the following way. In his second 
paper,! Samuel describes nitrous oxide as being 
formed from N2+O, that is, from an excited, 
divalent nitrogen molecule plus a normal oxygen 
atom with its two unpaired electrons. Since 
Samuel considers the excited nitrogen molecule to 
have the structure VIII, he is of course again led 
to the same classical structure V as before for the 
nitrous oxide. However, this structure does not 
have the same meaning in terms of eigenfunctions 
as it did before. If Samuel’s configuration of the 


| 
-N=N+— 
IX 


| 
N=N— 
VIII 


excited nitrogen molecule is examined in the 
manner employed above for that of the excited 
nitric oxide, it is found to correspond, from the 
modern point of view, to the structure 1X, rather 
than ‘VIII. Consequently, the resulting nitrous 
oxide must then be represented by the second 
modern structure VII. It is seen, therefore, that 
Samuel’s method of setting up an eigenfunction 
for the classical structure V can lead to either one 
of two quite different functions, depending upon 
whether the molecule is regarded as formed from 
NO+WN or from Ne+O. Moreover, these two 
different functions are just the ones which would 
have been set up to describe the two resonating 
modern structures VI and VII. 

We come now to consider the sense in which 
Samuel’s eigenfunction of excited nitric oxide, for 
example, can be said to correspond to structure 
III (as maintained by Samuel) rather than IV (as 
would follow from the modern theory). The point 
is that Samuel proposes a different way of 
correlating eigenfunction and structure. In par- 
ticular, he denies the existence of antibonding 
orbitals (cf. section C, below) and maintains 
instead that any orbital is bonding to a greater or 
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less degree if it is occupied by a pair of electrons 
contributed by different atoms. In this way, 
although he apparently retains most of the essen- 
tial features of the modern viewpoint, he uses.a 
different procedure for counting the number of 
bonds, and so he arrives at different conclusions 
in regard to structure. 

I shall not discuss in detail whether Samuel’s 
position in this respect is, or is not, a tenable one, 
but I should like to point out that, in practice, it 
is not always unambiguous. For example, there 
are two equivalent ways of regarding a system 
which consists of two helium atoms. From the 
pair-bond point of view, we would assign two 
electrons to each atom; since each electron is 
therefore localized about a single atom, and no 
electrons are shared, Samuel’s definition leads to 
the result that no bond exists. On the other hand, 
from the molecular-orbital point of view, we 
would assign two electrons to the lower (sym- 
metric) molecular orbital and two to the higher 
(antisymmetric) molecular orbital; since each 
orbital contains electrons contributed by differ- 
ent atoms, each is bonding, according to Samuel, 
and so a double bond exists. The present 
discrepancy is not due to any inconsistency 
between the pair-bond and molecular-orbital 
treatments, because the two corresponding eigen- 
functions prove to be identical. Consequently, if 
Samuel’s definitions are adopted, the eigenfunc- 
tion in question can be said to correspond to 
either zero or two bonds, depending entirely upon 
our mental attitude toward it. On the other hand, 
if we recognize that only the lower molecular 
orbital is bonding, whereas the higher is anti- 
bonding, then the discrepancy vanishes, since ro 
bond exists from either point of view. 

Let us now consider Samuel’s configuration of 
the excited nitric oxide molecule. As we saw 
above, it can be written in the form 


2507, 2sn, 2pmo', o?(2p), 2pmn?. 


Of all the orbitals involved, only the one desig- 
nated as o(2) contains electrons contributed by 
different atoms. Consequently, even when Samuel’ s 
own method of counting is used, there is only a 
single bond between the nitrogen and oxygen 
atoms, and so the structure of the excited nitric 
oxide molecule is really IV, rather than III as 
assumed by Samuel. It follows therefore that no 


way has as yet been defined for setting up an 
approximate eigenfunction to represent a struc- 
ture like III, and hence of a classical one like I 
or V, without the use of M-shell atomic orbitals. 

In concluding this section, I should like to call 
attention to the fact that Samuel’s arguments 
against the modern structures, which are based 
upon dipole moments, ionization potentials, and 
so on are irrelevant. No such experimental data 
can distinguish between two structures, one 
classical and one modern, which are merely two 
different ways of describing graphically the same 
eigenfunction. In particular, it may be noted that 
Samuel’s eigenfunctions for the classical struc- 
tures do not involve the configuration sp‘ of the 
neutral nitrogen atom as he states, but that 
instead, just like those for the modern structures, 
they involve the configuration sp* of the positive 
ion. 


C 


As was mentioned in section B, Samuel denies 
the existence of antibonding orbitals. To support 
this view, he brings forth spectroscopic evidence 
to show that a premoted® orbital is sometimes 
actually more strongly bonding than an un- 
premoted one. This argument is based upon the 
supposition that, in the molecular orbital theory 
of Mulliken et al., all premoted orbitals are con- 
sidered to be antibonding and all unpremoted 
orbitals to be bonding. As far as I know, however, 
no such view is held by any worker in the field at 
the present time.” The distinction between anti- 
bonding and bonding orbitals is based rather 
upon the presence or absence of nodes in the 
region between the two nuclei involved.® For 
example, let a and b be two atomic orbitals on 
different atoms, and let their signs be determined 
so that each is positive in the region where they 
overlap most strongly. Then, of the two molecu- 
lar orbitals, y and ¢, that are approximately 
representable as a+b and a—8, respectively, the 
former has one less node between the nuclei than 
the latter. It is considered, then, that y is 


6 R. S. Mulliken, Phys. Rev. 41, 49 (1932). 

7 See, for example, G. Herzberg, Molecular Spectra and 
Molecular Structure (Prentice-Hall, Inc., New York, 1939), 
page 398. ; 

8’ F, Hund, Zeits. f. Physik 73, 1 (1931); R. S. Mulliken, 
J. Chem. Phys. 1, 492 (1933). 
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ordinarily bonding whereas @¢ is ordinarily 
antibonding. 

The reasoning which leads to this supposition 
is as follows. An electron in the orbital y spends 
a larger part of its time in the region between the 
two nuclei than would one in either of the atomic 
orbitals a and 6, whereas an electron in the 
orbital ¢ spends a smaller part of its time there. 
Consequently, the average force of electrostatic 
attraction between the electron and the two 
nuclei, which is the only force capable of leading 
to the formation of a bond, must be greater in y, 
but smaller in ¢, than that in either a or 0b. In 
other words, an electron in y makes the approach 
of the two nuclei easier, but one in ¢ makes it 
more difficult, so that y is bonding whereas ¢ is 
antibonding. 

The foregoing argument indicates that the 
existence of both bonding and antibonding 
orbitals is a necessary consequence of the re- 
spective distributions of electric charge. It also is 
of interest in connection with Samuel’s further 
contention that the molecular-orbital treatment 
contains no theory of valence at all inasmuch as 
(according to Samuel) it ignores the factor re- 
sponsible for the bonding forces, namely, the 
exchange degeneracy of a pair of electrons. It 
may be noted also that, even in the pair-bond 
treatment, in which the exchange degeneracy is 
considered, and which Samuel therefore regards 
as the only actual theory of valence, the distri- 
bution of charge is still the decisive factor. In a 
system consisting of two helium atoms, for ex- 
ample, exchange degeneracy still exists, but there 
is no bond-formation because the charge-distri- 
bution is such that the atoms actually repel each 
other strongly—more strongly, in fact, than if 
there were no exchange degeneracy. It appears, 
therefore, that the exchange degeneracy is im- 
portant in the theory of valence only insofar as it 
affects the distribution of charge. If, as in the 
molecular-orbital treatment, a _charge-distri- 
bution leading toe attraction between the atoms 
can be arrived at without explicit consideration 
of the exchange degeneracy, the treatment loses 
none of its theoretical significance merely because 
it ignores this (for its purposes) irrelevant ex- 
change degeneracy. 
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D 


In this concluding section, I wish to mention 
two rather less fundamental matters which, how- 
ever, underly much of Samuel’s discussion of 
dipole moments. 

The first of these questions is concerned with 
Samuel’s calculation of the percent ionic charac- 
ters of essentially covalent bonds® between dis- 
similar atoms. He treats hydrogen chloride, for 
example, in the following manner. The dipole 
moment of this molecule can be considered due 
to a displacement of the two electrons that pro- 
duce the bond between the hydrogen and chlorine 
atoms. From the molecular orbital point of view, 
these electrons are in the same orbital, and so 
they must be displaced together. Now, the dis- 
placement of two electrons by the internuclear 
distance of 1.28A, as in Samuel’s calculation, 
leads to a moment of 24.8107" e.s.u. 1.28 
X 10-* cm = 12.2 D. Since the observed moment is 
only 1.03 D,* the percent ionic character is there- 
fore (1.03/12.2) X100=8.5 percent. 

Let us consider the significance of this method 
of calculation employed by Samuel. The structure 
I may be written to represent an imaginary state 
of the hydrogen chloride molecule, in which the 
bond is a pure covalent one with zero dipole 
moment. The average position of the two valence 
electrons is therefore halfway between those of 
the nuclei, and so 0.64A from either. Structure 


H—Cl H+Cl- H-Cl* 
I I] III 


II, on the other hand, may be written to repre- 
sent an equally imaginary state of the molecule, 
in which the bond is a pure ionic one, with a 
moment of 4.8107! e.s.u.X1.28X10-* cm 
= 6.1 D. Similarly, structure III may represent a 
second pure ionic structure, with a moment of 
6.1 D in a direction opposite to that of structure 
II. Now, the displacement of two electrons by the 
whole distance, giving rise to a moment of 12.2 D, 
would lead to the structure II only if one starts 
with III, and not if one starts with |. In order to 
go from II to I, one must instead displace the 
pair of electrons by only half the internuclear 


® Pages 182 and 183 of reference 1. 

* Actually, Samuel uses an incorrect value of 1.28 D for 
this moment, so that he ends up with about 10 percent ionic 
character, instead of 8.5 percent. 
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distance, that is, by 0.64A. The resulting moment 
is then, of course, again only 6.1 D, instead of 
12.2 D. Samuel’s method of calculation makes it 
appear as if he is under the obviously erroneous 
impression that III, rather than I, represents the 
pure covalent bond with zero dipole moment. 
The only way in which Samuel’s calculation 
can be given a self-consistent interpretation 
seems to be the following. If the three structures 
I, Il, and III make contributions of X percent, 
Y percent, and Z percent, respectively, to the 
dipole moment, and if W is the mean of Y and Z, 
then Y=W+8.5 and Z=W-—8.5. The figure 8.5 
is then taken as the percent ionic character. It 
would seem both simpler and more logical to 
follow Pauling,'® and to consider, as an extremely 
rough first approximation, that resonance with 
III may be ignored (i.e., that Z=0), and that, 
furthermore, the percent ionic character is given 
simply by Y, which is then equal to 8.5+8.5 
=17.0. It should be apparent that this same 
value of 17 percent ionic character is obtained, 
whether one employs the molecular-orbital treat- 
ment, as Samuel does, and considers the displace- 
ment of both electrons by only half the inter- 
nuclear distance, or employs the valence-bond 
treatment, as Pauling does, and considers the 
displacement of only one electron by the whole 
distance. It follows, therefore, that the dis- 
crepancy between Samuel’s and Pauling’s values, 
which Samuel considers a very important one 
resulting from a fundamental difference between 
the molecular-orbital and pair-bond treatments, 
is due merely to a difference in definition. 
Samuel’s percent ionic character is defined in such 


“L. Pauling, The Nature of the Chemical Bond (Cornell 
University Press, Ithaca, New York, 1939), first edition, 
p. 46. 


a way that it is always exactly half as great as 
Pauling’s. 

The second point to be discussed here is con- 
cerned with the fact that the moment of the 
N#1—0O bond, which is very nearly equal to zero 
in nitric oxide, is equal instead to about 2.0 D in 
a nitroso compound. Samuel’s explanation" of 
the variation is that the nitrogen atom has a 
single unpaired electron in nitric oxide, which it 
shares with a carbon atom, and which it therefore 
loses control of to some extent, in the nitroso 
compound. Consequently, a positive charge is 
produced on the nitrogen, and a dipole moment 
results. This explanation, which is used without 
essential change with several further types of 
molecule, ignores the following facts. First, the 
nitrogen atom, when it loses to the carbon atom a 
share of the electron which it formerly possessed, 
also gains from the carbon atom a share of an 
electron which the latter formerly possessed. If 
the sharing of these two electrons is equal, no 
dipole is produced. Second, any moment arising 
in the above manner would belong to the carbon- 
nitrogen, and not to the nitrogen-oxygen, bond. 
And finally, if, as a result of an inequality in 
sharing, the nitrogen atom does gain a positive 
charge, the resulting moment would be in the 
wrong direction; Samuel fails to realize this fact 
because he attributes the moment to the wrong 
bond. 
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On the Possibility of Estimating the Average 
Surface Area of Solid Particles 


Hans M. CAssEL 
General Ceramics and Steatite Corporation, Keasbey, New Jersey 
April 14, 1945 


T would be easy to derive the average surface area of 

crystalline powders from the thermal effects caused by 
submersion into a suitable liquid if the heat of emersion 
per unit area were known. Since this is not the case it 
appears as a very ingenious idea to coat the crystalline 
particles with the adsorbed film of a liquid whose surface 
energy can be independently determined. 

No objection shall here be raised against the principle of 
this method, but, unfortunately, the experimental proce- 
dure described by Harkins and Jura! implies a source of 
error which invalidates the reported results. 

With the objective of surface conditioning single crystal- 
lites, the contiguous mass of the powder is exposed to the 
saturated vapor of a liquid and the entire system (pre- 
sumably) allowed to attain thermal equilibrium. The 
powder treated in this way is then plunged into a body of 
the same liquid and the thermal effect observed. 

In evaluating the measurements, it is assumed that the 
contact angle at the liquid-crystal-vapor boundary is zero. 
However, the point has been overlooked that under this 
very condition capillary condensation between the con- 
tracting crystals is inevitable. Consequently, in a thus 
treated powder, groups of particles may stick together, 
held by the capillary forces of the liquid. However, even 
the immediate cohesion of adjacent crystals may cause 
inseparable clusters to form. 

Effects of this kind have been repeatedly investigated. 
Tomlinson? measured the adhesion between quartz fila- 
ments and concluded that the magnitude of the forces 
appears altogether too large for an explanation depending 
on the surface tension -of a liquid film. Stone* found that 
glass beads would not adhere unless the humidity is close 
to saturation, but once sticking together cannot be sepa- 
rated by passing dry air to remove the moisture. In the 
work of Bloomquist and Shutt‘* it has been shown that 
microscopic glass spheres suspended in organic liquids are 
flocculated on increasing the amount of dissolved water. 
The coalescence is correlated to the magnitude of the 
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interfacial tension at the organic liquid-water film bound- 
ary. Bradley’ demonstrated a striking increase in the sur- 
face tension of concentrated suspensions caused by the 
mutual attraction of neighboring solid particles. 

Taking all this into consideration, in the described pro- 
cedure an indefinite reduction of the immersed film area 
has to be expected. 

On the other hand, if capillary condensation is not ac- 
complished the assumption of perfect wettability cannot 
be true. In that case, the heat of emersion per unit area 
would be smaller than anticipated and, again, the surface 
area would be underestimated. 

Of course, such troubles would be avoided if it were 
possible to measure the thermal effects of particles indi- 
vidually coated and individually dropped into a liquid. 


1 William D. Harkins and George Jura, J. Am. Chem. Soc. 66, 1362 
) 


(1944). 
2G. A. Tomlinson, Phil. "1 6, 695 (1928). 
3 ‘tag Stone, Phil. Mag. 9, 610 (1930). 
4C. R. Bloomquist and R. S. Shutt, Ind. a Chem. 32, 827 (1940). 
5S. R. Bradley, J. Phys. Chem. 43, 663 (1939). 





Ion-Migration Phenomena Observed with the 
Electron Microscope 


HELMUT THIELSCH 
University of Michigan, Ann Arbor, Michigan 
April 30, 1945 


DROPLET of water containing dissolved KCI (5 g 

per liter) was placed on a collodion film and then in- 
cluded into a vacuum system which was pumped down to 
10-* mm Hg within 30 seconds. Viewing the result in the 
electron microscope three rather curious types of structures 
are apparent, a regular cubic crystal structure (Fig. 1), 
a regular linear dendritic system (Fig. 1), and an irregular 
branch structure (Fig. 2). Obviously, the differences are 
due to the rate of nucleus formation and growth, de- 
pending, of course, upon the variation of local ion con- 
centration within the droplet of water during the process 
of evaporation. 

The first mentioned type of crystal (the cubes shown in 
Fig. 1) results from a relatively slow rate of growth; the 
second linear type will be produced by a soméwhat faster 
rate (its closeness to the more stable cubic crystals signi- 
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fies that), and the irregular dendritic branch structure re- 
quires a still faster rate. 

The two dendritic structures are unstable in atmospheric 
conditions, since upon standing for 6 hours in the atmos- 
phere (Figs. 3 and 4), or merely breathing upon the 
sample (Figs. 5 and 6) will cause migration. In the last 
two figures outlines of the previous structure are still 
visible. 


Whether or not the streaks or schlieren which appear 
in regularly grown crystals (Fig. 7) are identical with the 
linear dendrites observed or are due to electric surface 
charges,* further investigations will have to bring out. 

A detailed discussion of this will follow in the near future. 


*H. Seemann, “Die wandernden Schlieren in elektronenoptischen 
Bildern von Einkristallen, insbesondere Molybdanoxyd,”’ Naturwiss. 31, 
415-416 (1943). 
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The Effect of Temperature on 
Molecular Parachors* 


Tuomas S, GARDNER 
Tennessee Eastman Corporation, Kingsport, Tennessee 
April 26, 1945 


HE effect of temperature on molecular parachors was 

not included in Professor Samuel’s discussion. The 

temperature effect is of sufficient magnitude to invalidate 

the use of molecular parachors for structural determina- 
tions for many organic compounds.! 

In order to use molecular parachor constants for struc- 
tural work, it would be necessary to standardize all data 
on a constant temperature basis. It was found! in the case 
of the organic compounds investigated that nearly all of 
them exhibited an increase in molecular parachor with 
temperature. This was interpreted as being due to dissocia- 
tion with increasing temperature. In the case of ortho, 
meta, and para isomers of compounds able to form chelated 
rings in the ortho form, the meta and para isomers gave the 
familiar effect of increasing parachor with increasing tem- 
perature, while the ortho isomer yielded almost constant 
molecular parachors. 

We are in substantial agreement with Professor Samuel 
on the impossibility of assigning a parachor value to the 
coordinate covalent double bond. 

The difficulty of determining parachor values for struc- 
tural work is illustrated in Table I by the comparison 


TABLE I 








Samuel’s Values 
(P) (P 


Buehler, Gardner, et al@ 
(P) 


Observed 
T(corr.) (P) 


30.0 
40.0 
50.1 


-Calcu- 
lated¢ 


Calcu- Ob- 
lated? served 


195.3 





196.8 196.2 
196.7 


198.6 


Sulfuryl Chloride 


40.0 
70.3 
100.9 
130.4 


297.3 
299.7 
301.2 
301.5 


Ethyl] Sulfite 


473.9¢ 
473.1 
473.8 


Diphenyl Sulfone 


189.6 
189.5 
190.6 
190.8 


Ethyl! Nitrate 


371.0 
372.0 
373.2 
375.5 


p-Toluene Sulfony! 370.0 


Chloride 








@ See reference 1. 

> Calculated by Samuel using parachor values for different valences. 

© Calculated using Sugden’s values for the parachor constants except 
that Sugden’s value of —1.6 for the coordinate bond has not been 
included. 

4 Samuel's calculated Sugden values allowed Sugden’s value for the 
coordinate bond. 

* Results for diphenyl sulfone were unpublished* because of the 
unexplainable, high anomaly. The analysis of the material was: melting 
point (literature) 128-129°C; found 127.5-127.9°C. Percent S: calcu- 
lated 14.70; found 14.90 and 14.94. The constancy of the parachor 
values for different temperatures is interesting. Unpublished data used 
in determining the observed parachors are as follows: 


T(corr.) D Y 
150.4 1.1434 38.08 
174.8 1.1256 35.54 
191.6 1.1139 34.27 
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of a few values obtained by both investigators. 

It is evident that Professor Samuel’s uses of parachor 
values for different valence states cannot compensate for 
the temperature effects of the parachors of associated 
organic compounds. 

* An addendum to Professor R. Samuel's discussion (J. Chem. Phys. 
12, 167 (1944)). 


1C, A. Buehler, T. S. Gardner, and M. L. Clemens, J. Org. Chem. 2, 
167 (1937). 





Comments on “The Effect of Temperature on 
Molecular Parachors’’* 


R. SAMUEL 
154 Nassau Street, New York, New York 
May 10, 1945 


HE parachor is indeed not entirely independent of 

temperature. Other molecular constants which are 
less sensitive to geometrical factors behave, however, in a 
similar manner, and this effect therefore cannot be ascribed 
entirely to association. For instance, the increased popula- 
tion of higher vibrational levels will probably be found 
among its causes. The various modes of vibration will be 
differently affected according to their different energies, 
which may result in a deformation of the molecule and may 
influence the intermolecular potentials. 

However, for reasons which are discussed in a forth- 
coming publication, no rigorous additivity can exist for 
any of the so-called additive molecular constants, and an 
additivity in first-approximation is all what can be ex- 
pected. Dr. Gardner’s figures show that with decreasing 
temperature the parachors indeed converge towards those 
values for which this first-order additivity has been es- 
tablished. It is very interesting to see that such values 
are already obtained at temperatures midway between 
melting and boiling points. 

The atomic parachors of course have been calculated for 
such values only, and these are, therefore, the only ones 
which can be used for structural work. 


* Preceding letter by Thomas S. Gardner. 





Remark to Dr. Wheland’s Article 


R. SAMUEL 
154 Nassau Street, New York, New York 
May 19, 1945 


HE article of Dr. Wheland appearing in this issue of 

this journal cannot be answered without detailed 
considerations of experimental evidence, the meaning of 
additivity, and so forth. As such a detailed answer could 
not be prepared in time for this issue for outside reasons, 
it will be found in a forthcoming paper. 

In this moment attention should be drawn only to the 
main point, namely, Dr. Wheland’s theoretical argument. 
The wave function he uses is that of a diatomic and not 
that of a polayatomic molecule with its increased number 
of orbitals. It has no bearing on the problem of pentacova- 
lent nitrogen. 





